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This paper gives details of the analytical and numerical procedures used 
to solve the basic problem of the scattering of a plane electromagnetic wave 
by an axisymmetric raindrop. A nonperturbative solution is obtained by 
expanding the scattered and transmitted fields in terms of spherical vector 
wave functions, so that Maxwell's equations are satisfied exactly in the 
regions exterior and interior to the raindrop, and by combining point 
matching with least-squares fitting to satisfy the boundary conditions on 
the surface of the raindrop with sufficient accuracy. 

Numerical results are presented for scattering by oblate spheroidal rain- 
drops, with eccentricity depending on (and increasing with) drop size, 
for two orthogonal polarizations of the incident wave. The calculations 
were made at 4, 11, 18.1, and 30 GHz, in the case in which the direction of 
propagation of the incident wave is perpendicular to the axis of symmetry 
of the raindrop, which is of interest for terrestrial microwave relay systems. 
At 30 GHz, the calculations were also made for the case in which the angle 
between the direction of propagation and the axis of symmetry is 70° and 
50°, since different elevation angles are of interest for satellite systems. 
These basic results were summed earlier over the drop-size distribution to 
calculate the differential attenuation and differential phase shift caused 
by rain, which are of importance in the investigation of cross polarization 
in radio communication systems. 
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We also derive the first-order perturbation approximation to the scatter- 
ing by axisymmetric raindrops that are nearly spherical, which generalizes 
Oguchi's results for spheroidal raindrops icith small eccentricity. Some 
simplifications that may be made in his formulas are pointed out. The 
perturbation results serve as a useful check on the least-squares-fitting 
procedure applied to spheroidal raindrops with small eccentricity. In 
addition, considerable improvement is obtained in the closeness of the 
perturbation results to the least-squares-fitting ones, in particular for the 
larger drop sizes, by perturbing about an equivolumic spherical raindrop, 
with appropriate perturbation parameter, rather than perturbing about an 
inscribed spherical raindrop, as did Oguchi. Similar comparisons were 
also made earlier for the rain-induced differential attenuation and differ- 
ential phase shift, and these quantities were calculated approximately at 
frequencies up to 100 GHz, using the results corresponding to perturbation 
about the equivolumic spherical raindrop. The perturbation results are 
obtained quite inexpensively, whereas the least-squares-fitting procedure is 
very costly. 

I. INTRODUCTION 

In a recent short note, the authors and Chu 1 gave calculated results 
of differential attenuation and differential phase shift caused by rain, 
based on scattering of a plane electromagnetic wave by oblate sphe- 
roidal raindrops. These results are of importance in the investigation of 
cross polarization in radio communications systems. In this paper, we 
give details of the analytical and numerical procedures used to solve 
the problem of scattering by a single raindrop, which were only out- 
lined in the note. Although the results given in this paper are for 
oblate spheroidal raindrops, the procedures are applicable for axisym- 
metric raindrops that are not too nonspherical, and calculations could 
be made for raindrops that are more flattened on the bottom than on 
the top, such as for the shapes determined by Pruppacher and Pitter. 2 

Two polarizations of the incident wave, designated I and II, are 
considered, as depicted in Fig. 1. The factor *"*•" has been suppressed. 
In the first polarization, the electric field is parallel to the plane con- 
taining the axis of symmetry of the raindrop and the direction of 
propagation of the incident wave. In the second polarization, the elec- 
tric field is perpendicular to this plane. The angle between the direction 
of propagation and the axis of symmetry is denoted by a. In terrestrial 
microwave relay systems, a = 90° is of main concern, but other values 
of a are of interest for satellite systems. 
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The incident wave induces a transmitted field in the interior of the 
raindrop and a scattered field. In the far field, the quantities of pri- 
mary interest are the complex forward scattering amplitudes, 3 <S'i(0) 
and (Sn(O). In the two polarizations considered, the polarization of the 
far scattered field in the forward direction is the same as that of the 
incident wave. Also of basic interest are the cross sections of the rain- 
drop. The total cross sections Q] and Q} 1 are given in terms of the real 
parts of the forward scattering amplitudes by eq. (39), where ko is the 
free space wave number. We also calculate the scattering cross sections 
Q] and Q] 1 . The absorption cross sections Ql and Q" are given in terms 
of the total and scattering cross sections by (38). 

In Section II, we discuss the formulation of the problem of the 
scattering of a plane electromagnetic wave by a single raindrop. 
Spherical coordinates are chosen with polar axis along the axis of sym- 
metry of the raindrop, and origin interior to it, as in Fig. 2. The scat- 
tered field is expanded in terms of solutions of the vector wave equa- 
tion, with wave number ko, which satisfy the radiation condition. An 
analogous expansion is assumed for the transmitted field inside the 
raindrop, in terms of vector wave functions, with wave number ki 
= Nk , which are finite at the origin. Here, N is the complex refrac- 
tive index of the raindrop. The complex coefficients in the expansions 
are to be determined by satisfying the boundary conditions, namely, 
the continuity of the tangential components of the total electric and 
magnetic fields across the surface of the raindrop. 

In Section III, the incident field is expanded in a (complex) Fourier 
series in the azimuthal angle <p. Because of the axial symmetry of the 
raindrop, the problem can be decomposed and the boundary conditions 
satisfied independently for each term of the Fourier series. In Section 
III, expressions are also given for the forward scattering amplitudes 
and the scattering cross sections in terms of the coefficients in the ex- 
pansion of the scattered field. In addition, we express the total and 
scattering cross sections for an elliptically polarized incident wave in 
terms of those for the two linearly polarized incident waves under 
consideration. 

In Section IV, we give expressions for the first-order approximations 
to the coefficients in the expansions of the scattered and transmitted 
fields for axisymmetric raindrops that are nearly spherical. These 
results generalize those given bj- Oguchi 4 for spheroidal raindrops with 
small eccentricity. Since our derivation follows closely the one given 
by Oguchi, we omit most details. However, we point out some simplifi- 
cations that may be made in his expressions. 
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In Section V, we discuss an approximate nonperturbative solution 
to the problem of scattering by an axisymmetric raindrop. For each 
term in the Fourier series expansion in the azimuthal angle <p, the four 
boundary conditions should be satisfied on the cross-sectional bound- 
ary curve r = R(8), ^ ^ ir, which defines the shape of the rain- 
drop. Only a finite number of coefficients in the expansions of the 
scattered and transmitted fields is considered. These coefficients are 
determined approximately by requiring the boundary conditions to 
be satisfied in a least-squares sense at a total number of points on the 
cross-sectional curve that is greater than the number of unknown 
coefficients. 

In Section V, we also discuss the advantage of using least-squares 
fitting rather than collocation, in which the total number of fitting 
points is equal to the number of unknown coefficients that are then 
determined by solving a system of simultaneous linear equations. 
After we had completed the calculations for scattering by oblate 
spheroidal raindrops at 4, 18.1, and 30 GHz, a paper was published by 
Oguchi 5 in which he carried out similar calculations for a = 90° at 
19.3 and 34.8 GHz and used collocation for the expansions in terms of 
spherical vector wave functions. At 34.8 GHz, he also used an expansion 
in terms of spheroidal wave functions and truncated the infinite system 
of equations which he derived from the boundary conditions. 

In Section VI, the least-squares-fitting program and some subsidiary 
programs are discussed. The numerical routines used for calculating the 
special functions that enter into the boundary conditions are also 
described. In addition, some indication of the running times involved 
and the storage requirements are given. 

In Section VII, we first discuss the checks that were made on the 
least-squares-fitting program. These include comparison with the Mie 
theory 6 of the results of scattering by spherical raindrops at different 
angles of incidence. We also compare extrapolated results for oblate 
spheroidal raindrops with small eccentricity to the first-order perturba- 
tion results. 

We then discuss our calculations of the scattering by oblate sphe- 
roidal raindrops, for which the ratio of minor to major semiaxis depends 
linearly on the radius a (in centimeters) of the equivolumic spherical 
drop ; specifically, 

a/b = (1 - a), ab 2 - a 3 . (1) 

This relationship is similar to that used by Oguchi. 4,5 The rain-induced 
attenuation and phase shift were calculated 1 for both polarizations by 
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summing the real and imaginary parts of the forward scattering ampli- 
tudes over the Laws and Parsons drop-size distribution, as quoted by 
Setzer. 7 Thus, for rain rates up to 150 milimeters per hour, there are 
14 different drop sizes, a = 0.025(0.025)0.35, to be considered. 

The calculations were done for wavelengths of 7.5, 2.727, 1.6575, 
and 1.0 cm, corresponding to frequencies of approximately 4, 11, 18.1, 
and 30 GHz. The refractive indices N at 20 °C were obtained from an 
elaborate fitting equation in a recently published survey 8 of available 
measured data (except at 4 GHz, for which older data were used, since 
the calculations at that frequency were made at an earlier date). The 
angle of incidence a was taken to be 90° at 4, 11, and 18.1 GHz, while 
at 30 GHz the calculations were done for a = 70° and a = 50° also. 
The calculated values of the forward scattering amplitudes Si (0) and 
*Sii (0) are given in Tables II to VII, and those of the total cross sections 
Q\ and Q} 1 and the scattering cross sections Q] and Q" are given in 
Tables VIII to XIII. Section VII concludes by discussing some cal- 
culations using collocation and mentioning some checks on Oguchi's 
calculations at 19.3 and 34.8 GHz. 

In Section VIII, we compare three sets of first-order perturbation 
results with those obtained by least-squares fitting for oblate sphe- 
roidal raindrops. One set of results corresponds to perturbation about a 
spherical raindrop with radius equal to the length a of the minor semi- 
axis of the spheroidal raindrop, which was the procedure used by 
Oguchi. 4 The other two sets correspond to perturbations about the 
equivolumic spherical raindrop, with different perturbation param- 
eters that are consistent for small drop sizes. Considerable improve- 
ment is obtained in the closeness of the perturbation results to the 
least-squares-fit results by perturbing about the equivolumic spherical 
raindrop with the appropriate perturbation parameter. The compari- 
sons are presented graphically in Figs. 3 to 14. The values of the 
forward scattering amplitude S(0) and the total and scattering cross 
sections Q, and Q, for the equivolumic spherical raindrops are given in 
Tables XIV to XVII. These values are independent of the polarization 
of the incident wave and of the angle of incidence a. 

The three sets of perturbation results for the differential quantities 
Q" - Q\ and Im[5i(0) - &n(0)] are compared in Figs. 15 to 23 with 
those obtained by least-squares fitting for oblate spheroidal raindrops. 
Again, considerable improvement is obtained by perturbing about the 
equivolumic spherical raindrop. In a recent short note, 9 similar com- 
parisons were made for the rain-induced differential attenuation and 
differential phase shift. Moreover, these quantities were calculated 
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approximately at frequencies up to 100 GHz, using the two sets of 
results corresponding to perturbations about the equivolumic spherical 
raindrop. The perturbation results are obtained quite inexpensively, 
whereas the least-squares-fitting procedure is very costly. 

The appendices contain some details that it was considered desirable 
to omit from the main text. 

II. FORMULATION OF PROBLEM 

We now consider the problem of scattering of a plane electromag- 
netic wave by a single raindrop. Suppressing the factor e~ iul , where w 
is the angular frequency, the divergenceless electric and magnetic 
fields E and H satisfy Maxwell's equations 10 



V X E = t'wMoH, V X H = (a - ta»e)E, 



(2) 



where mo is the constant permeability, a is the conductivity, and e is 
the dielectric constant. Exterior to the raindrop a = and e = e , 
while interior to it a = <j\ and e = ci. The appropriate boundary con- 
ditions 11 are that the tangential components of the total electric and 
magnetic fields be continuous across the surface of the raindrop. Let 



k 2 = a>/to(a>e + ia), 



(3) 



with Re(fc) > 0. Then the free space wave number is k = cof>oeo)* 
and the wave number in the raindrop is 



ki = Nk , 



(4) 



where N is the complex index of refraction of water. 

We consider two polarizations of the incident wave depicted in Fig. 
1 . We choose Cartesian coordinates (x, y, z) with origin interior to the 



EjjeikoUsina * zcosal 



S n <0)E n - 




Fig. 1 — Two polarizations of the incident wave. 
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raindrop and z-axis coinciding with the axis of symmetry of the rain- 
drop. The direction of propagation of the incident wave is perpendicu- 
lar to the y-axis and inclined at an angle a to the z-axis. In the first 
polarization, the magnetic field is assumed parallel to the y-axis and 
the incident fields are given by 

Ei = £i(cos ai — sin ak) exp [ik (x sin a + z cos a) J, 

fco " ' (5) 

Hi = — Eij exp [ik (x sin a + z cos a)l. 

cono 

where i, j, and k denote unit vectors parallel to the coordinate axes. In 
the second polarization, the electric field is assumed parallel to the 
y-&xis and the incident fields are given by 

En = En j exp [iko(x sin a + z cos a)] 
and (6) 

— kn 

Hit = .#11 (cos ai — sin ak) exp [ik (x sin a + z cos a)]. 

We now consider the problem of representing the scattered and 
transmitted fields induced by the incident wave. It is convenient to 
introduce spherical coordinates (r, 6, tp) with corresponding unit vec- 
tors ii, i 2 , and i 3 as depicted in Fig. 2. Then the equations 

V X M = m, V X N = fcM (7) 

are satisfied by the spherical vector wave functions, 12 

M..(*) = z, (*r)«<- '[ & Pi"' (cos tiU ~ " P l "'j;° S A U] (8) 
and 
*».(*) = g—L(n + 1) !^ Pi" 1 (cos e)u + [ 2 -^ + *(*!■)] 

X [^W) ia + jL^co S ,),]j. (9, 

Here z„ denotes a spherical Bessel function 13 of order n and Pj, m| de- 
notes the associated Legendre function 14 (of the first kind) of degree n 
and order \m\, where m is a positive or negative integer, and n is an 
integer with n i: \m\ and n ^ 0. The prime denotes derivative with 
respect to the argument. As a matter of convenience, we have chosen 
to use complex linear combinations of the even and odd spherical 
vector wave functions. 12 
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Fig. 2 — Cartesian and spherical coordinates. 



Outside the raindrop, the total electromagnetic field is the sum of 
the incident field of the plane wave and the scattered field. The scat- 
tered field must satisfy the radiation condition and, consequently, in 
view of eqs. (2), (3), and (7), we assume expansions of the form 



E>= - Z E [a»,Mffi(* ) + 6»nNffl(*o)] (10) 

m = — » ng \m\ 

715*0 



and 



H- = *• £ E La m Jf&(ko) + 6..M®(fc )l (11) 

"MO m = -« nS|m| 

where the superscript 3 denotes that spherical Bessel functions of the 
third kind, i.e., spherical Hankel functions of the first kind, are used. 
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Thus, in (8) and (9), z„(k r) = h ( n l) (k r). For k r » 1, 

^ 1) (M~^r^ 1 fi aor , (12) 

so that the expansions in (10) and (11) involve outgoing waves. 

Analogous expansions are assumed for the transmitted field inside 
the raindrop except that, since the origin of the coordinate system is 
interior to the raindrop, spherical Bessel functions of the first kind 
must be used so that the field remains finite at r = 0. Also, the wave 
number inside the raindrop is k\, as given by (4). Thus, we assume ex- 
pansions of the form 



E' = - L £ [c»„M^(fci) + rf mn N^(A;i)] (13) 



m = — »o ?i2|m| 
n ^0 



and 



H' = **i ± Z Cc-irffiObi) + <*™Mffl(*i)l (14) 

WMO m=- =o n^\m\ 



where the superscript 1 indicates that z n (kir) = j n (kir) in (8) and (9). 
The unknown (complex) coefficients a mn , b m „, c m „, and d mn in (10), 
(11), (13), and (14) must be determined from the boundary conditions. 
The surface of the raindrop is given by 

r = R(6), ^ e ^ v, ^ <p ^ 2w, (15) 

where it is assumed that R(6) is a single- valued, continuously differ- 
entiable function of 8. The continuity of the tangential components of 
the total electric and magnetic fields across the surface of the raindrop 
implies that, for r = R(6), 

El + Ei = E' 3 , (16) 

m + m - m, (17) 

« + « + !§« + « -« + |g*, (18) 
771 + m + hw m + H ° affi + sf H{ > (19) 

where 2£y = Eiy and Hj = Hiy and the incident fields E* and H* 
are given by (5) or (6) . 
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III. FAR-FIELD QUANTITIES 

Because of the axial symmetry of the raindrop, it is convenient to 
expand the incident plane wave in a (complex) Fourier series in the 
azimuthal angle <p, and we write 

E*= f e«(r, 0)e**", H* = £ h m (r, B)e im *. (20) 



It follows from (5) and (6) that 



k 



and 



ei(r, 6) = EJ m (r, 6), K>(r, 6) = E l -i g m (r, 5) (21) 

CO/io 

eg(r, 0) = tfng»(r, 5), hS(r, 0) = -tfn -^ f-(r, 0), (22) 



where expressions for f m (r, 0) and g„,(r, 0) are derived in Appendix A 
and are given by eqs. (79) and (80). 

If the boundary conditions (16) to (19) are multiplied by er im * and 
integrated with respect to <p from to 2ir, then a set of four equations 
involving the unknown coefficients a mn , b m „, c m „, and d mn is obtained 
for each m. These equations are given by (81) to (84) in Appendix A, 
where we have used the notations e m ; = e m ij and h mj = h m -iy. It 
foUows readily from (21), (22), and (79) to (84) that, for the first 
polarization of the incident wave, 

tf = M 

(23) 



0-— mn O'mny 


h 1 = h 1 
<J—mn — "msj 


Cf—mn = Cmri) 


ft 1 — H 1 


polarization 




a" = a 11 

** — mn u mni 


h u — —ft 11 


r II _ r U 
<s—mn t/mn 


"■—mn "■mn- 



(24) 

Thus, it is sufficient to consider only nonnegative values of m . 

It is worth noting that, if the raindrop is symmetrical about the 
plane 6 = tt/2, i.e., R (it - d) = R(0),0 ^ 6 ^ tt/2, then some further 
reductions may be made. In particular, in the case a = ir/2, it is found 
that 

dm,\m\+2s+l = = C m ,\m\+2s+lt . . 

|T _ n _ 71 \^) 

O m ,\m\+2s — U — a m> | m |+2si 

and 

0>m,lm\+2s = = C m ,|m|+2si . . 

&m,|m|+2s+l = = dm,\m\+2s+l, 

for s = 0, 1, 2, • • •, so that alternate coefficients vanish. Reductions 
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may still be made in the case ofa^ ir/2 by considering the sum and 
the difference of the boundary conditions corresponding to a and to 
<$ = (ir — a), as shown in Appendix B. We have not utilized these 
reductions in the program for calculating the unknown coefficients, 
since we wanted the program to be applicable to raindrops without a 
meridional plane of symmetry, that is, those raindrops flattened more 
on the bottom than on the top. However, (25) and (26) served as a 
useful check on the program in the calculations for spheroidal rain- 
drops with a = ir/ 2. 

We describe in Section V how we obtain approximate values of (a 
finite number of) the coefficients a mn , b m „, c mn , and d mn , but we now 
turn to the quantities of physical interest. We consider only the far 
scattered field, so that fc r»l. Thus, we restrict our attention to the 
leading term in the asymptotic expansion of the spherical Bessel func- 
tion of the third kind, as given by (12). Also, it follows that 

A " 1),(M ~ ^T eik ° T - (27) 

Then, from (8) to (11), it is found that 

knT6~ '* or E* 

~ t Z (-0-{a-.[^^l,-^rt-'(cos»)i,l 

m=-=on§|m| ( . «# Sin J 

~ *»■[ "^JF* «• + £i ^ (™ <*] ) «*" < 28 > 

and 

comoH 8 ~ kdi X E*. (29) 

Of particular interest are the scattered fields in the forward direction, 
corresponding to = a, <p = 0, for which, from (76), 

(cos oi — sin ak) =12, j = U- (30) 

From (5), (6), (23), (24), and (28) to (30), it follows that the far 
scattered field in the forward direction has the same polarization as the 
incident wave for either polarization. The forward scattering ampli- 
tudes are 3 

fljfl)) = i- (cosai - sin ok)- lim { -ik re- ikor El\e^ a , „_oi (31) 



and 



«M0) = J- j- lim {—ik<pe- ik °<E\i\$-.. v =o}. (32) 

•Oil r— =° 
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Thus, for the first polarization of the incident wave, 

EiSi(0) = £ L (-0- 1 

m = — oo n2|m| 

x[ai n .^Pl»'(co3a) + M. dPl" 1 («"■»> 1 (33) 
L sin a aa I 

and for the second polarization, 

^ii^ii(o) = £ e (-t)* +2 

m = — oo ;, g m 

X a%n -y + b% n • -. PJT 1 (cos a) . (34) 

|_ da sin a J 

The energy scattered by the raindrop is 15 

W. = ^ Re r [' VBi(HO* - Ei(Hi)*y sin 9ddd<p, (35) 

^ Jo Jo 

where the asterisk denotes complex conjugate. The calculation of W„ 
using the asymptotic form of the scattered fields given by (28) and 
(29) and letting r — > <x> t is outlined in Appendix C. It is found that 

W. = \ t E "(" + ))("+ j m |)j (| „,|. + | 4 .,|.). (36) 

a>MO&Om=T±=o n i\m\ (2tt + 1) (W - |m|)I 

The scattering cross section Q, is defined as the ratio of the scattered 
energy flow to the mean energy flow of the incident wave per unit area. 
Thus, 16 

r _ jWg Q u _ 2"MoWT ,a 7) 

W koEiSi' ^ ~ koEnETn {6i) 

The total (extinction) cross section is the sum of the scattering and 
absorption cross sections, so that 

Q] = Q\ + Ql, QV = QV + Ql 1 . (38) 

(We note that van de Hulst 16 uses the notations C ezt , C, ca , and C„6. for 
Qi, Q„ and Q a , respectively.) It is "known" 17 that 

Q\ = % ReSi (0), Q? = fj ReSn (0), (39) 

so that (38) may be used to determine the absorption cross sections Ql 
and Ql 1 . The relations (39) which are consistent with the optical 
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theorem may be verified directly from the relations 

Ql _ 2-ntin a n _ jjgWjg . . 

and the expression for the total energy 15 

W, = 1 Re /J' ^ [^(#2)* + #!(#')* 

- Si(ff|)* - El{Hi)*y* sin Od»d^, (41) 

and a few of the details are given in Appendix C. 

Let us now consider an (in general) elliptically polarized incident 
wave that is the sum of the two linearly polarized incident waves 
under consideration; i.e., 

E' = Ei' + Ei, (42) 

Then the scattered electric field is 

E* = Ef + Eix, (43) 

and, as shown in Appendix C, 

q _ (EiEiQ] + EnEliQl 1 ) .... 

(EjEj + EuEn) 



and 



Qt = (giggj + E u E*nQV) , (45) 

(EiEi + EuEii) 



Since for polarizations I and II the far scattered field in the forward 
direction has the same polarization as the incident wave, it follows from 
(31), (32), and (43) that the far scattered field in the forward direction 
for the elliptically polarized incident wave is given by 

E'|«-«. , =0 ~~^[tfi£i(0)(cosai - sinak) + £ n S„ (0)j]. (46) 

Thus, from (39) and (44) to (46), it suflices to calculate £i(0), &n(0), 
Qs, and Q n . The relation between the polarizations of the incident field 
and the far scattered field in the forward direction may be determined 
from (42) and (46)., using (5) and (6). 

IV. FIRST-ORDER PERTURBATION THEORY 

Oguchi considered spheroidal raindrops with small eccentricity and 
carried out a perturbation expansion originally determining the first- 
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order approximation 4 and later the second-order one. 18 We have cal- 
culated the first-order approximation for axisymmetric raindrops that 
are nearly spherical, so that the surface of the raindrop is given by 
r = R(0),0 ^ 6 ^ x, where 

R(fi) = <z[l + valid) +•••], l"l « 1- (47) 

Our derivation follows closely the one given by Oguchi. 4 Since the cal- 
culations are somewhat lengthy and involved, we merely outline the 
procedure, state the results, and point out some simplifications that 
may be made in the expressions given by Oguchi. 

The incident wave may be expanded in terms of spherical vector 
wave functions, 419 and the expansions are given by eqs. (116) and (117) 
in Appendix D. These expansions are consistent with those given by 
Oguchi, but the reader should bear in mind that, in addition to using 
the even and odd spherical vector wave functions, Oguchi has as- 
sumed the time factor e +,ut , and his waves propagate in the opposite 
direction to ours. Corresponding to (47), the coefficients in the expan- 
sions (10), (11), (13), and (14) are expanded in the form 

a mn = a£ + va™ + • • •, b mn = ftfi + 4fi + • ■ -, (48) 

c mn = d°> + wcSL + • • -, d mn = <&°> + vd™ + • • -. (49) 

Appendix D indicates how these coefficients may be determined from 
the boundary conditions (81) to (84). 

The zero-order approximation, with v = in (47), corresponds to a 
spherical raindrop of radius a. We have, for n ^ \m\ and n^O, 

offi = a mn a n , b™ = ro „&„ (50) 



and 
where 



Cmn — OtmnCti) (*mn — PmnCtnt v."-*-/ 



-i^(2n+D(n- |m|)l m fl-i (oOB tt) _ «. f (52 ) 



Ex n(n+ \){n+ \m\)\ sine* v B u 

PL = -i» +1 (2n + l)(n - \m\)\ dPl ml (cos a) = ^"_„ 
Ei n(n + l)(n + \m\)\ da En 



(53) 



and expressions for the quantities a n , b n , c„, and d n (which do not depend 
on the polarization) are given by eqs. (119) to (122) in Appendix D, 
where p = k a and the functions F„(£) and £„(£) are defined in (118). 
For a = 0, the coefficients vanish unless \m\ = 1, and the well-known 
Mie solution 6 is recovered. 
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The first-order corrections to the coefficients are given by 

a& = j n (Np)X mn , c™ = h™(p)X mn , (54) 

WS = G n (Np)Y m „ + j n (N P )Z mn , (55) 



and 



where 



dill = F n ( P )Y mn + I *?(p)Z Mi (56) 



X mn = (1 - JV*)p»c E Ld^Gi{Np)JJH + icWMNpVZl, (57) 

I £ I m I 

F m „ = (JV 2 - l) P 3 d n E M ji(iW3 - id%lGi(Np)I$2, (58) 

and 

Z m „ = -z(iV 2 - l)pd„ E dSHMNp)HS. (59) 

The quantities H n n h /$, and J,™ in (57) to (59) involve integrals over 
the perturbation of the raindrop surface from the sphere. Specifically, 

2(n + | w | ) ! „ m _ 



(2n + l)(n - |iw|)! 

= Z(Z + 1) /"' Pj m| (cos 0)Pl m| (cos 0) sin da 1 (d)dd, (60) 
y o 

2n(n + l)(^+ |wi|)l Jm _ ,„ 
(2n + l)(n- \m\)\ * nl 

[* [ rfPj m| (cosfl) rfPj, 7 " 1 (cos g) 

~Jo L do dd 

+ -^T« P ^ m| ( cos e)P ' m| (cos ^ I sin ^W*' ( 61 ) 

hill 1/ | 



and 

2n(w + l)(n + Iml)! 



« nl — c)nl 



(2n + l)(n - |m|)! 

dPJ, m| (cos 0) 



m m r rpj m| (cose) 



rfPJ m| (cosfl) 
dd 



Pi m| (cos0) \n(e)de. (62) 



The above results were obtained after considerable algebra and 
after using the simplifications given in eqs. (125) to (131) in Appendix 
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D. The expressions in (54) to (56) hold for both polarizations of the 
incident wave, and the only quantities therein that depend on the 
polarization are the zero-order coefficients c$' and dj$, which enter into 
the expressions given in (57) to (59) and are given by (51) to (53). 
In general, there are infinitely many terms in the sums in (57) to (59) , 
but in particular cases there are only a finite number of terms. 
Thus, for a spheroidal raindrop with 

R(0) = o(l - v sin 2 0)-l = o(l + \v sin 2 d + ■ • •), (63) 

we have <ri(d) = \ sin 2 0, from (47). The integrals in (60) to (62) may 
be evaluated explicitly, and it is found that HZ and 1% vanish unless 
I = n, n + 2, or n — 2, and J$ vanishes unless I = n — 1 or n + 1. The 
explicit expressions for these quantities are given by eqs. (137), (139), 
and (142) in Appendix E, in which 8 in denotes the Kronecker delta, 
i.e., 8in = 1 for I = n, and otherwise. We have verified that our results 
for the spheroid are consistent with those of Oguchi, 4 provided that 
simplifications corresponding to those in eqs. (123) and (127) to (130) 
are made in his expressions, and due allowance is made for the differ- 
ences in notation. We remark that similar simplifications may be made 
in Oguchi's expressions even in the case in which the permeability of the 
spheroid differs from the free space value. 

As is seen later, it is advantageous to obtain the first-order approxi- 
mation for a spheroidal raindrop by perturbing about the equivolumic 
sphere, rather than the inscribed sphere as Oguchi 4 did. If a is the 
radius (in centimeters) of the equivolumic sphere, then, from (1), 

a = a(l - a)*, v » o(2 - a). (64) 

Hence, from (63), 

R(8) = a[_l + 2a(£ sin 2 B - \) + 0(a 2 )] 

= a[l + K!sin 2 0-i)+O(* 2 )]. (65) 

We must now replace p by p = k d and add terms to the expressions in 
(60) to (62) corresponding to 5i(d) = — f . It is readily found, using the 
orthogonality relation for the Legendre functions 14 and equations (102) 
and (103), that these additional terms correspond to 

8% = -\n(n + 1)5,», 1% = -\B ln , 7S = 0. (66) 

We also remark that the use of the perturbation parameter v = 2a, 
rather than u as given by (64), generally gives better results for the 
larger drop sizes. 

970 THE BELL SYSTEM TECHNICAL JOURNAL, JULY-AUGUST 1 974 



V. LEAST-SQUARES-FITTING PROCEDURE 

We now discuss the calculation of an approximate nonperturbative 
solution of the scattering problem. As mentioned in Section III, it is 
sufficient to determine the unknown coefficients for nonnegative values 
of m and then to use the relationships in (23) or (24). For m = 0, 1, 2, 
• • • , the boundary conditions in (81) to (84) take the form 

K mq (e) - T, \_a mn A mnq {9) + b mn B mnq {9) 

n £m 

+ c mn C mng (d) + d mn D mnq (6)l = 0. (67) 
for q = 1, 2, 3, 4 and ^ 9 ^ t, where 

K ml (9) = e m3 (R(9), 9), K m2 (9) = ^-° h m3 (R(9), 9), (68) 
K m3 (9) = e m2 (R(9), 9) + ^ ^ e ml (R (6), 9), (69) 

K mi (9) = ^ \h mi (B(e),e) +]^^*»itfW)»')} ( 7 °) 



and 



The function R(9) describes the shape of the raindrop. The functions 
e mj = e m -i y and h mj = h m i y are given by (21) or (22), depending on 
the polarization of the incident wave, where expressions for f m (r, 9) 
and gjf(r, 6) are given by (79) and (80). The functions A mnq {9), 
B mnq (6), C mnq (9), D mnq {9), which do not depend on the polarization, 
involve the spherical Bessel functions of the first and third kinds and 
the associated Legendre functions and the derivatives of each of these 
functions. In view of (4), the argument of the spherical Bessel functions 
of the first kind is complex. 

For each m there are infinitely many unknown coefficients a mn , b mn , 
c m „, and d mn . To obtain an approximate solution, only a finite number 
of coefficients is considered. One procedure is to truncate the sum in 
(67) atn = No, say, and then to satisfy the boundary conditions at 
the points 9 = 9 tm , 1=1, ■•-, (AT — m + 1 - 5 m0 ), which are ap- 
propriately selected, e.g., uniformly spaced in the interval to tt. This 
was the procedure adopted by Oguchi, 5 and it leads to a system of 
simultaneous linear equations for the coefficients. We refer to this pro- 
cedure, in which the total number of fitting points is equal to the 
number of unknown coefficients, as collocation. 

The method of collocation was used by Mullin et al. 20 for the much 
simpler two-dimensional scalar problem of scattering by a perfectly 
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conducting cylinder of smooth contour, which is not a gross perturba- 
tion from the circular. In this problem, there is only one set of unknown 
coefficients to be determined, namely, that occurring in the expansion 
of the scattered held. Mullin et al. checked the results for the circular 
cylinder obtained by collocation with the known analytical results and 
also considered elliptic cylinders (with ratio of minor to major axis of 
|, in particular). Before tackling the raindrop problem, we considered 
the same problems as Mullin et al., but combined point matching with 
least-squares fitting. 

Thu.s, instead of using collocation, we satisfied the (single) boundary 
condition in a least-squares sense at a larger number of points than the 
number of unknown coefficients in the truncated expansion of the 
scattered field. We found that a significant improvement could be 
obtained in the overall fit of the boundary condition, although the far 
field quantities were not affected as significantly. This is because the 
higher-order coefficients are more significant on the boundary than in 
the far field. However, the accuracy of the lower-order coefficients is 
affected by the goodness of fit of the boundary condition. With colloca- 
tion, there were much larger errors in the boundary condition (in 
between the fitting points) than with least-squares fitting with a suffi- 
ciently large number of points. 

Since the fit of the boundary condition for the elliptic cylinder be- 
comes poorer with increasing eccentricity, we considered it desirable 
to use least-squares fitting rather than collocation for the raindrop 
problem. Thus, in order to approximately satisfy the boundary condi- 
tions (67), we minimized the quantities 

A m = £ w mq £ \K mq (e, m ) - £ [a m »A nnq {0i m ) 

9=1 1=1 n=m 

+ b mn B m n„(e lm ) + c m „c mnq (e lm ) + d mn D mnq (d lm )2\ 2 (71) 

for each m = 0, •• •, M, with respect to the (complex) coefficients a mn , 
b mn , c m „, and d mn , where w mq > are appropriate weights and 0* m are 
appropriate points in the interval to ?r. It is assumed that 

X m ^ N m - m + 1 - 8 m0l (72) 

so that the total number of fitting points is not less than the number 
of unknown coefficients to be determined. In the case of equality in 
(72), least-squares fitting is equivalent to collocation. 

The programs for carrying out least-squares fitting and for calculat- 
ing the special functions occurring in the functions in (71) with argu- 
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merit &i m are described in the next section. Actually, more flexibility was 
built into the least-squares-fitting program, allowing for truncation of 
the sums in (67) at different limits for each of the coefficients a m „, b mn , 
c mn , and d mn , and for X m and 6i m in (71) to depend on q, so that each of 
the four boundary conditions could be fit at different points, and in 
particular at a different number of them. It was anticipated that the 
least-squares-fit subroutine might become overloaded, in which case 
it would be desirable to hold the number of coefficients to a minimum. 
It turned out, however, that the subroutine was able to handle almost 
100 (complex) coefficients without difficulty. Similarly, it might be 
desirable to keep the total number of fitting points to a minimum, and 
hence to use fewer fitting points for those of the four boundary condi- 
tions that are easier to fit. Again, it was not found necessary to do this 
for the calculations carried out so far. 

For the calculations of this paper we took the weights to be inde- 
pendent of m and q, i.e., w mq = 1, since it was generally found that the 
difference was tolerable between the magnitudes of the maximum error 
in the fit of each of the four boundary conditions. We at first considered 
factoring out (sin 0) m_l , for m ^ 2, from the boundary conditions (67), 
but decided against it since we felt that the absolute, rather than the 
relative, error in the fit of the boundary conditions was important. 
However, because of the presence of this factor, we did experiment 
with unequally spaced fitting points which were closer together in the 
neighborhood of ir/2. We decided that equally spaced fitting points 
would suffice, provided that enough were used. The total number of 
fitting points was usually taken to be slightly more than twice the num- 
ber of unknown coefficients, i.e., X m > 2(N m — m -f- 1 — 8 m0 ). 

Generally, N m , the upper limit of n in (71), was taken to be inde- 
pendent of m, i.e., N m = N , m =0, ■ ■-, M. The choice of N and M 
depended both on drop size and on frequency. The choice of M was 
based on the rate of convergence of the outer sums in the expressions 
in (33), (34), and (36) for the far field quantities. On the other hand, 
to ensure the accuracy of the lower-order terms in the inner sums, it is 
necessary to take more terms in n than are really needed in the cal- 
culation of the far field quantities. A convergence test was carried out 
by doing the calculations for N m = No, (N + 2) and (N + 4). 

VI. NUMERICAL ROUTINES 

The program to compute the complex coefficients a mn , b m „, c mn , d mn , 
the scattering cross section Q a , and the forward scattering amplitude 
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S(0) for the two polarizations of the incident wave is written in Fortran 
IV for a Honeywell 6070 computer. It uses the complex arithmetic and 
math routines (such as sin, cos) written for that system. The program 
is written as a three-part package; the first part consists of a driver 
routine which sets up core storage for the least-squares matrix and 
associated vectors and the following subroutines : 

L2FIT — The main subroutine which computes the elements of the 
complex matrix for input to the least-squares-fitting pro- 
cedure and controls the other subroutines. 
CLSTSQ — Computes a least-squares fit for a linear system with com- 
plex coefficients; the algorithm and Fortran routine were 
written by P. Businger of Bell Laboratories. 
BJYNC — Computes the Bessel functions J n (z), Y n (z), for z complex, 
n a nonnegative integer; the algorithm and Fortran sub- 
routine package were written by E. Sonnenblick of Bell 
Laboratories. 
SBES — Computes the spherical Bessel functions j n (x), y„ (x), x real, 
n a nonnegative integer. 
CSB — Computes the spherical Bessel function j n (z), z complex, 
n a nonnegative integer. 
SLEG — Computes the associated Legendre functions P™(cos0), m, 
n nonnegative integers. 

The second part of the package is the routine 

SQS — Computes the forward scattering amplitude *S(0) and the 
scattering and total cross sections Q s , Qt for both polariza- 
tions from the least-squares-fit solutions. 

The third part is a computational check on the least-squares fit. It 
consists of a driver routine as in the first part, the function subroutines 
BJYNC, SBES, CBS, SLEG, and the main subroutine 

CHECK — Checks the accuracy of the least-squares computation of 
a>mn,b mn , c mn , d mn by using these coefficients to compute the 
boundary fit at points in addition to those used to obtain 
the coefficients. 

The internal computations of all the subroutines except CSB (see 
detailed description below) are done in double-precision arithmetic; 
on the Honeywell 6070 this means 18 digits are used for all computa- 
tions. The function values (Bessel, spherical Bessel, Legendre), how- 
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ever, are returned as single-precision numbers (8 digits), since the 
accuracy of any more digits could not always be guaranteed. 

The main subroutine, L2FIT, sets up the complex matrix to mini- 
mize the quantities (71). Instead of the limit N m in (71), the program 
actually breaks the summation on re into four parts, using limits a m , 
m , 7m, 5 m for the coefficients a m „, b mn , c mn , d mn , respectively; further, 
it replaces the limit X m in the sum on I with four limits, \ mq . Thus, for 
each m = 0, 1, • • •, M the matrix equation to be minimized takes the 
form ||A-x — b|| £^ 0, where A is an L-by-iV" matrix with L, N, which 
are both dependent on m, defined as follows: 

4 

Li = £_, "mq 

3=1 

and 

N = a m + /3 m + y m + 8 m - 4(??* - 1 + 5 m0 ). 

Although the basic functions (Bessel, Legendre, and derivatives of 
these) are returned to L2FIT in single precision, the calculations of the 
elements A mng , B mnqt C mng , D mnq in the least-squares matrix and K mq 
in the vector of constants, b, are carried out and left in double precision 
for input to the least-squares-minimization program. To facilitate 
changing the raindrop shape and spacing of points on the boundary, 
the quantities d tm and R(d [m ) are computed in subroutines called by 
this routine. 

The routine CLSTSQ uses elementary Hermitian (or Householder) 
transformations to compute a linear least-squares solution to the 
equation ||Ax — b|| = min. The algorithm is an adaptation of an 
algorithm for a real matrix written by P. Businger and G. H. Golub. 21 

The routine CHECK uses the coefficients a m „, b m „, c mn , d mn from 
CLSTSQ as input into the expressions in (67) to compute the fit of 
the boundary condition in between the fitting points as well as at the 
fitting points. The goodness of this fit provides a check on the accuracy 
of the computed coefficients. 

The computation of the elements for the least-squares matrix and 
constant vector requires values of J n (x), j n (x), y n (x), j n (z), PJ?(cos 9), 
where x and 6 are real, z is complex, and n and m are nonnegative 
integers. The routine to compute J„(x), BJYNC 22 uses a downward 
recursion scheme (compute J N , Jn-\, ■ ■ •, Jo) for |x| > 0.1 and uses 
the power series expansion of J„(x) for \x\ ^ 0.1. Comparison of the 
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results of this routine with tables in Abramowitz and Stegun 23 for x 
values ranging from to 5 yielded complete agreement in all cases. 

The real spherical Bessel functions, j n (x) and y n (x), are computed 
recursively in the routine SBES using the algorithm of D. S. Drum- 
heller, 24 an improved version of Miller's algorithm. 

Miller's algorithm for j n (x) is: Let g N +\(x) = 0, g N (x) = 10" 20 (some 
small fixed number); generate g N -i(x), •', Go(x) from the recurrence 
relations, g n -i(x) = (2n + l)xr l g n (x) - g n +i(x), satisfied by the 
spherical Bessel functions; 25 compute j (x) = xr 1 sin x ; normalize ji(x) 

Drumheller's algorithm generates ascending orders of y n (x) re- 
cursively, starting from Rayleigh's formulas 26 for y (x), yi(x) up to 
some order N, where N is strictly limited only by the relation 

y N (x) ^ max — the largest number the computer can handle 
(for the Honeywell 6070, max « 10 38 ). 

Letting 

In(x) m 0, / M (») = -{xh, N {x))-\ 

the algorithm then generates f N -%{%), ' ' ', Mx) recursively, using the 
recurrence relations satisfied by the spherical Bessel functions. Al- 
though, as Drumheller points out, f (x) = jo(x) to some degree of 
precision (determined to a large extent by N), to ensure exactness in 
the lower orders of j n (x) and to shift any error to the higher-order, 
smaller-magnitude terms, we calculate j (x) exactly and normalize 
ji(x), • • •, j/r(x) as in Miller's algorithm. We compared these results 
to tables in Abramowitz and Stegun 26 and in the U. S. Math Tables 
Project 27 for the values of n and the range of x used in the least-squares 
fitting, and the agreement was excellent. 

The complex spherical Bessel function, j n (z), is computed in the 
routine CSB, using an algorithm designed by A. E. Kaplan, 28 since 
Drumheller's (or Miller's) recursive algorithm yields inaccurate results 
for complex arguments with a significant imaginary part. Kaplan uses 
a Taylor series expansion to compute jw(z), Jn-i(z), with N tt |z| 2 
for best results, then uses the backward recursion and normalization 
techniques discussed above to generate /nr-i(f)i ■ ■ ', /•(*)■ Tne use of 
the Taylor series for "large" complex arguments produces better start- 
ing values and, therefore, more accurate recursion results than either 
Drumheller's or Miller's algorithm. A "large" complex argument in 
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this case is one such that \z\ ^ 9 ; for \z\ < 9, we use Drumheller's 
algorithm directly. In this routine, complex arithmetic is used for all 
internal computations ; thus, since the Honeywell 0070 does not have 
double-precision complex arithmetic, all computations are in single 
precision. The accuracy of the results of this routine was checked by 
the following comparisons : 

(i) For z real, they were compared to the spherical Bessel functions 

in Abramowitz and Stegun 26 and to results from SBES. 
(ii) For z pure imaginary, they were compared (multiplied by ap- 
propriate factors) to the modified spherical Bessel functions 
of the first kind in Abramowitz and Stegun. 29 
(in) For z any complex number, they were compared to Rayleigh's 
formulas 26 of order 0, 1, 2 (a straight computation of Rayleigh's 
formulas produces inaccurate values for higher orders). 

For all values of z (from | z | = to | z | = 50) and for order n as large 

as we could compute (or as large as we could compare to in tables or 

formulas), this routine produced answers agreeing in six to eight 

decimal places with the other results. 

The associated Legendre function P£(cos0) is computed directly 

from its series expansion, in powers of sines or cosines, as given by 

L. Robin. 30 It should be noted here that we use the following definition 

for PJ? (cos 0), 

d m P n (cos 6) 



P£(cos0) = (-l) m sin m 



d(cosd)' 



whereas Stratton 31 omits the factor ( — l) m . Comparison of the results 
of this routine to explicit formulas for m ^4, n ^ 7 yielded eight 
places of agreement. Further checks of this routine against tables given 
by S. L. Belousov 32 and tables in the U. S. Math Tables Project 33 were 
done by S. Hoffberg for the values of m and n used in the least-squares 
fitting ; her results also found complete agreement in all cases. 

In the scattering problem, the matrix and vector sizes, as well as 
the highest order needed for the functions, depend on drop size and 
frequency ; the sizes required by these parameters are discussed in the 
next section. Here, we give approximate running times for each of the 
three parts of the package and give some feeling for the correlation of 
the least-squares matrix size and the limit max m with the overall core 
storage and run time in the first part. In (71), we generally took N m 
m No, and X m dependent upon N and decreasing with m (for m ^ 1). 
Then the largest matrix size occurs when m = 0, and is thus a function 
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of No alone. Tabulated below are some typical storage and run time 
figures for the least-squares fit. 



max m 


No 


6 

7 
8 


13, 15, 17 
15, 17, 19 
19, 21, 23 



Largest 
Matrix Size 

164 X 68 
180 X 76 
212 X 92 



Approximate Approximate 
Storage Run Time (hours) 



70 K 

80 K 

103 K 



0.13 

0.2 

0.39 



In each of the above cases, we computed a set of coefficients a mn , b m „, 
c nn , d mn , n = m, ■ •-, No (n 9* 0), form - 0, 1, •••, max m, and for 
each of the three values of No. 

The second and third parts require much less storage because neither 
requires a large storage matrix. The total core requirement to run SQS 
is 12 K; it computes the quantities S(0), Q„ Qt in typically less than 1 
second. The third part uses 35 K to do the calculations for 428 rows and 
92 columns (more rows are needed because the boundary fit is checked 
at 8 values in addition to those used to obtain the coefficients) ; to 
check the L2FIT results for N = 23 required approximately 0.05 
hours. 

As a final note, we save the coefficients computed from L2FIT on 
magnetic tape. At first, all coefficients were written on tape and the 
second and third parts used the tape as input; this proved very in- 
efficient, due in part to the high cost of tape usage and in part to the 
fact that we were saving all data generated. We switched to writing 
the data from each run of L2FIT onto a high-speed disc file, using this 
as input to the other two parts ; this change resulted in a noticeable 
cost reduction and allowed us to permanently save only the best data. 



Table I — Raindrop parameters for different drop sizes 



a(cm) 


o(cm) 


V 


0.025 


0.02458158 


0.049375 


0.05 


0.04831913 


0.0975 


0.075 * 


0.07120149 


0.144375 


0.1 


0.09321698 


0.19 


0.125 


0.11435330 


0.234375 


0.15 


0.13459757 


0.2775 


0.175 


0.15393617 


0.319375 


0.2 


0.17235477 


0.36 


0.225 


0.18983822 


0.399375 


0.25 


0.20637045 


0.4375 


0.275 


0.22193444 


0.474375 


0.3 


0.23651206 


0.51 


0.325 


0.25008398 


0.544375 


0.35 


0.26262956 


0.5775 
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Table II — Forward scattering amplitudes at 4 GHz 
with a = 90° for different drop sizes 



o(cm) 



Si(0) 



Sn(O) 



0.025 

0.05 

0.075 

0.1 

0.125 

0.15 

0.175 

0.2 

0.225 

0.25 

0.275 

0.3 

0.325 

0.35 



6.9215 X 10~ 8 
5.8893 X 10' 7 
2.2523 X 10-o 
6.3705 X 10-" 
1.5414 X 10- 6 
3.3841 X 10- 6 
6.9433 X 10- & 
1.3555 X 10"« 
2.551 X 10-" 
4.68 X 10- 4 
8.45 X 10-" 
1.51 X lO- 3 
2.72 X 10- 3 
4.9 X 10- 3 



r« — 8.6909 X 10"H 



6.8473 X 10"H 
2.2822 X 10~H 
5.3582 X 10-H 
1.0399 X 10"H 
1.7918 X 10~H 
2.8484 X 10"H 
4.276 X 10'H 
6.154 X 10-»» 
8.59 X WrH 
1.170 X 10~H 

■ 1.57 X 10-H" 
2.08 X 10-H 

- 2.7 X 10-H 



7.3309 X 10- 8 
6.5803 X 10" 7 
2.6370 X 10- 6 
7.7684 X 10- 6 
1.9518 X 10~ 5 
4.4527 X 10- 6 
9.5324 X 10- 5 
1.9565 X 10- 4 
3.916 X 10- 4 
7.77 X 10~< 
1.553 X 10- 3 
3.20 X 10-» 
6.96 X 10"* 
1.63 X 10-* 



8.9487 X 10-H' 

- 7.2647 X 10-H' 
2.4970 X 10-H' 
6.0510 X 10-H 

- 1.2134 X 10-H 
2.1630 X 10-H 

- 3.5628 X 10-H 

- 5.552 X 10-H 

- 8.318 X 10-H 
- 1.212 X 10-H 
- 1.735 X 10-H 

- 2.46 X 10-H 

- 3.47 X 10-H 

- 4.8 X 10-H 



VII. LEAST-SQUARES-FITTING RESULTS 

We begin by discussing the different ways in which the least-squares- 
fitting program was checked. First, calculations were carried out at 
both 4 and 34.8 GHz for centered spherical raindrops, corresponding 
to R(6) = a in (15). The results were compared with the calculations 
based on the zero-order solution given in Section IV, corresponding to 
v = in (47). Comparison was made for several values of a and differ- 
ent values of a, and excellent agreement was obtained for the far-field 
quantities, generally to six or seven significant figures. As expected, 
the far-field quantities are independent of the angle of incidence a. 

Table III — Forward scattering amplitudes at 11 GHz 
with a = 90° for different drop sizes 



a (cm) 


S 


i(0) 


Sn(O) 


0.025 


4.8189 X 10" 6 


- 1.8194 X 10-H 


5.0841 X 10~ 6 - 


- 1.8734 X 10-H 


0.05 


5.9423 X 10" 6 


- 1.4657 X 10-H 


6.5120 X 10" 6 - 


- 1.5555 X 10-H 


0.075 


3.6151 X 10- 4 


- 5.0737 X 10-H 


4.0920 X 10-* - 


- 5.5577 X 10-H 


0.1 


1.6675 X 10- 3 


- 1.2532 X 10-H 


1.9652 X 10- 3 - 


- 1.4186 X 10-H 


0.125 


6.5377 X 10" 3 


— 2.5280 X 10-H' 


8.0766 X 10" 3 - 


- 2.9419 X 10-H 


0.15 


2.0109 X 10-2 


— 4.0292 X 10-H 


2.4674 X lO- 2 - 


- 4.6783 X 10-H 


0.175 


3.7203 X 10- 2 


- 4.7914 X 10-H 


4.1732 X 10-2 . 


- 5.8328 X 10-H 


0.2 


4.735 X 10-* 


- 5.734 X 10-H 


5.492 X 10-2 . 


- 8.098 X 10-H 


0.225 


5.958 X 10-* 


— 7.504 X 10-H 


7.929 X 10-2 - 


-1.1621 X 10-H 


0.25 


7.67 X 10" 2 


- 9.68 X 10-H 


1.173 X 10- 1 - 


- 1.580 X 10-H' 


0.275 


9.89 X lO- 2 


— 1.225 X 10-H 


1.725 X 10-' - 


- 2.060 X 10-H 


0.3 


1.29 X 10' 1 


— 1.51 X 10-H 


2.51 X 10" 1 - 


- 2.54 X 10-H 


0.325 


1.67 X 10-» 


— 1.78 X 10-H' 


3.54 X 10- 1 - 


- 2.93 X 10-H 


0.35 


2.1 X 10- 1 


— 2.0 X 10-H 


4.8 X 10" 1 


-3.1 X 10-H 
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Table IV — Forward scattering amplitudes at 18.1 
with a = 90° for different drop sizes 



GHz 



a (cm) 


-Si(O) 


Sn(O) 


0.025 


4.0158 X 10- 6 


— 8.1579 X 10"H 


4.2246 X 10' 5 — 8.4000 X 10-H 


0.05 


6.5425 X 10-" 


— 6.7265 X 10-H" 


7.1168 X 10~ 4 —7.1424 X 10"H 


0.075 


5.0674 X 10' 3 


— 2.3476 X 10-H 


5.6959 X 10-' — 2.5717 X 10~H 


0.1 


2.2608 X 10- 2 


— 5.1254 X 10~H 


2.5696 X 10~ 2 — 5.7588 X WrH 


0.125 


5.0374 X 10- 2 


— 8.0587 X 10-«i 


5.7722 X 10- 2 — 9.6353 X 10-H 


0.15 


8.5403 X 10~ 2 


— 1.2201 X 10~H' 


1.0834 X 10- 1 — 1.5714 X 10~H 


0.175 


1.3950 X lO" 1 


— 1.7392 X 10-H' 


1.9921 X 10-' — 2.3098 X 10~H' 


0.2 


2.1700 X 10-' 


— 2.2871 X 10-H 


3.3903 X 10-' — 2.9883 X 10"H 


0.225 


3.181 X 10- 1 


— 2.742 X 10"H 


5.229 X 10-» — 3.320 X 10~H 


0.25 


4.307 X 10-' 


— 2.999 X 10-H 


7.173 X 10- 1 —3.136 X WrH 


0.275 


5.390 X 10-» 


— 3.089 X 10-H' 


8.899 X 10"' — 2.633 X 10~H 


0.3 


6.36 X 10" 1 


— 3.12 X 10-H 


1.038 — 2.14 X 10-H" 


0.325 


7.22 X 10-' 


— 3.20 X 10-H 


1.179 — 1.84 X 10-H 


0.35 


8.0 X 10- 1 


— 3.4 X 10-H 


1.34 — 1.8 X lO'H 



Moreover, the values of the coefficients obtained from least-squares 
fitting were checked against those calculated from (50) and (51), 
subject to (52), (53), and (119) to (122). 

Next, the least-squares fit was carried out for spherical raindrops 
when the origin of the coordinate system was offset from the center of 
the raindrop, so that 



R(d) = aid cos 6 + (1 - 8 2 sin 2 0)*]. 



(73) 



The calculations were done for different values of a, 8, and a, with the 
largest value of 8 being 0.325 at 4 GHz and at 0.2 at 34.8 GHz. As ex- 



Table V — Forward scattering amplitudes at 30 GHz 
with a = 90° for different drop sizes 



a (cm) 


5i(0) 


S,i(0) 


0.025 


3.4513 X 10- 


* — 3.7481 X 10-H 


3.6235 X 10- 


* — 3.8595 X 10-H 


0.05 


6.9873 X 10- 


3 — 3.1187 X 10-H 


7.5859 X 10- 


3 — 3.3144 X 10-H 


0.075 


4.5783 X 10- 


2 — 9.5267 X 10-H 


5.1071 X 10- 


2 — 1.0490 X 10-H 


0.1 


1.3415 X 10- 


1 — 1.8677 X 10-H' 


1.6165 X 10- 


1 — 2.1613 X 10-H 


0.125 


2.9755 X 10- 


1 — 2.8388 X 10-H 


3.8979 X 10" 


1 — 3.2171 X 10-H 


0.15 


5.1727 X 10- 


1 — 3.3781 X 10-H 


6.9041 X 10- 


i — 3.3561 X 10-H 


0.175 


7.3731 X 10- 


1 — 3.4278 X 10-H 


9.6534 X 10" 


» — 2.7438 X 10-H 


0.2 


9.274 X 10"' 


— 3.461 X 10-H 


1.2001 


— 2.302 X 10-H' 


0.225 


1.1122 


— 3.885 X 10"H 


1.4661 


— 2.419 X 10-H 


0.25 


1.3309 


— 4.693 X 10-H 


1.8221 


— 2.627 X 10-H 


0.275 


1.601 


— 5.57 X 10-H 


2.245 


— 2.21 X 10-H 


0.3 


1.902 


— 6.23 X 10-H 


2.662 


— 1.20 X 10-H 


0.325 


2.20 


— 6.70 X 10-H 


3.06 


— 2.4 X 10-H 


0.35 


2.49 


— 7.3 X 10-H 


3.50 


+ 4 X 10-H 
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Table VI — Forward scattering amplitudes at 30 GHz 
with a = 70° for different drop sizes 



o(cm) 


S.(0) 


Sn(0) 


0.025 


3.4679 X 10- 


« — 3.7610 X 10~H 


3.6200 X 10- 4 


— 3.8594 X lO'H 


0.05 


7.0267 X 10- 


3 — 3.1417 X 10-H' 


7.5553 X 10- 3 


— 3.3145 X 10-H' 


0.075 


4.6185 X 10" 


2 — 9.6618 X 10-H 


5.0856 X 10- 5 


— 1.0513 X 10-4 


0.1 


1.3701 X 10- 


1 — 1.9131 X 10-H' 


1.6133 X 10" 


— 2.1727 X 10-H 


0.125 


3.0867 X 10" 


1 — 2.9154 X 10"H 


3.9039 X 10- 


— 3.2498 X 10-H 


0.15 


5.4041 X 10- 


1 — 3.4346 X 10-H 


6.9413 X 10- 


— 3.4113 X 10"H 


0.175 


7.7133 X 10" 


' — 3.4313 X 10-H 


9.7408 X 10" 


— 2.8097 X 10-H' 


0.2 


9.723 X 10-' 


— 3.425 X 10-H 


1.2143 


— 2.357 X 10-H' 


0.225 


1.1720 


— 3.832 X 10-H 


1.4838 


— 2.454 X 10-H' 


0.25 


1.4123 


— 4.621 X 10 _ H 


1.8412 


— 2.678 X 10-H 


0.275 


1.710 


— 5.45 X lO'H 


2.269 


— 2.34 X 10-H 


0.3 


2.042 


— 6.04 X 10-H" 


2.699 


— 1.41 X 10-H 


0.325 


2.38 


— 6.46 X 10-H 


3.11 


— 4.8 X 10-H 


0.35 


2.71 


— 7.0 X 10-H 


3.56 


+ 1.5 X 10-H 



pected, the far-field quantities are independent of 8, as well as a, and 
again excellent results were obtained. These calculations provided a 
nontrivial check on the programming of the boundary conditions, 
since dR/dd f£ 0. In addition, they gave some idea of the increase in 
the number of terms in n that is required, a result of the ratio of 
maximum to minimum distance from the raindrop surface to the 
origin, which is necessarily greater than unity for oblate spheroidal 
raindrops. 

As a final check on the least-squares-fitting program, calculations 
were carried out at 34.8 GHz for oblate spheroidal raindrops with small 

Table VII — Forward scattering amplitudes at 30 GHz 
with a = 50° for different drop sizes 



a(cm) 


Si(0) 


Sn(O) 


0.025 


3.5101 X 10- 


* — 3.7936 X 10-H 


3.6111 X 10" 


4 — 3.8590 X 10-H 


0.05 


7.1265 X 10" 


3 —3.1998 X 10-H 


7.4779 X 10- 


3 — 3.3147 X 10 H 


0.075 


4.7204 X 10- 


2 — 1.0005 X 10-H 


5.0312 X 10- 


2 — 1.0571 X 10-H 


0.1 


1.4431 X 10- 


1 — 2.0284 X 10"H 


1.6054 X 10" 


i — 2.2016 X 10-H' 


0.125 


3.3714 X 10- 


' — 3.1102 X 10-H 


3.9193 X 10- 


1 — 3.3329 X 10-H' 


0.15 


6.0012 X 10- 


1 — 3.5745 X 10" H 


7.0369 X 10" 


i — 3.5523 X 10-H 


0.175 


8.5967 X 10" 


i — 3.4211 X 10-H 


9.9679 X 10- 


1 — 2.9777 X 10-H 


0.2 


1.0888 


— 3.278 X 10-H 


1.2518 


— 2.491 X 10-H' 


0.225 


1.3247 


— 3.580 X 10-H' 


1.5306 


— 2.520 X 10-H 


0.25 


1.6145 


— 4.252 X 10-H 


1.8902 


— 2.752 X 10-H 


0.275 


1.974 


— 4.91 X 10-H 


2.327 


— 2.58 X 10-H 


0.3 


2.379 


— 5.27 X 10-H 


2.785 


— 1.87 X 10-H' 


0.325 


2.80 


— 5.45 X 10-H' 


3.24 


— 1.03 X 10-H 


0.35 


3.24 


— 5.7 X 10-H' 


3.72 


— 4 X 10-H 
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Table VIII — Total and scattering cross sections at 4 GHz 
with a = 90° for different drop sizes 



a (cm) 


Qftcm)' 


Q?(cm)* 


Q](cmy 


QViwx)* 


0.025 


1.2393 X 10"« 


1.3126 X 10" 6 


8.9950 X 10" 10 


9.5367 X 10" 10 


0.05 


1.0545 X 10- 6 


1.1782 X 10- 5 


5.5428 X 10-" 


6.2399 X lO" 8 


0.075 


4.0327 X 10~ 6 


4.7215 X 10- 5 


6.0811 X IO- 7 


7.2816 X lO" 7 


0.1 


1.1406 X 10-* 


1.3909 X lO' 4 


3.2922 X 10-6 


4.2012 X 10-" 


0.125 


2.7599 X 10' 4 


3.4947 X 10- 4 


1.2106 X 10-5 


1.6500 X 10- 8 


0.15 


6.0592 X 10- 4 


7.9725 X 10- 4 


3.4869 X 10" 6 


5.0884 X 10"' 


0.175 


1.2432 X 10-' 


1.7068 X 10~ 3 


8.4904 X 10" 6 


1.3303 X lO" 4 


0.2 


2.427 X 10- 3 


3.503 X 10" 3 


1.830 X 10-* 


3.088 X 10- 4 


0.225 


4.568 X 10- 3 


7.011 X 10- 3 


3.599 X 10" 4 


6.567 X 10" 4 


0.25 


8.38 X 10-" 


1.391 X lO" 2 


6.60 X lO-* 


1.309 X lO" 3 


0.275 


1.51 X 10" 2 


2.78 X lO" 2 


1.15 X 10-* 


2.50 X 10- 3 


0.3 


2.71 X 10" 2 


5.73 X 10- 2 


1.93 X lO' 3 


4.7 X 10- J 


0.325 


4.87 X 10" 2 


1.246 X 10-> 


3.2 X lO" 3 


8.8 X lO" 3 


0.35 


8.8 X lO" 2 


2.92 X 10" 1 


5.2 X lO" 3 


1.78 X lO' 2 



eccentricity, corresponding to v = 0, 0.05, 0.1, and 0.15 in (63). The 
calculations were done for a = 90° and for a = 0.025(0.025)0.275. 
Corresponding to (48), the total cross section may be expanded in the 
form Q t = Ql 0) + vQl l) + ■ ■ -. Values of Q} a) and Qj ul) were obtained 
from the least-squares results by extrapolation and were compared with 
the perturbation values given by Oguchi. 4 Unfortunately, there were 
significant discrepancies for the larger drop sizes, the largest error 
being more than 17 percent. Consequently, we did the perturbation 
calculations ourselves and obtained results differing from our extrap- 



Table IX — Total and scattering cross sections at 11 GHz 
with a = 90° for different drop sizes 



o(cm) 


Q](cmy 


& n (cm) 2 


Qj(cm) 2 


Q s "(cm) 2 


0.025 


1.1407 X lO" 5 


1.2035 X lO" 5 


5.1556 X lO" 8 


5.4666 X lO" 8 


0.05 


1.4066 X 10" 4 


1.5415 X lO' 4 


3.2095 X lO- 8 


3.6172 X lO" 6 


0.075 


8.5573 X 10" 4 


9.6863 X 10" 4 


3.5933 X 10-* 


4.3203 X lO" 5 


0.1 


3.9471 X 10~ 3 


4.6519 X lO" 3 


2.0259 X 10- 4 


2.6114 X lO" 4 


0.125 


1.5476 X lO" 2 


1.9118 X 10" 2 


8.1236 X 10" 4 


1.1302 X lO" 3 


0.15 


4.7600 X lO" 2 


5.8407 X lO" 2 


2.7238 X lO" 3 


4.0701 X lO" 3 


0.175 


8.8064 X 10- 2 


9.8785 X lO" 2 


7.3803 X 10- 3 


1.1535 X lO" 2 


0.2 


1.1208 X 10" 1 


1.3001 X lO" 1 


1.559 X 10- 2 


2.643 X 10" 2 


0.225 


1.4103 X 10- 1 


1.8769 X lO -1 


2.815 X 10- 2 


5.325 X 10- 2 


0.25 


1.815 X 10" 1 


2.777 X 10- 1 


4.60 X lO" 2 


9.76 X 10- 2 


0.275 


2.341 X 10"' 


4.084 X 10" 1 


7.21 X 10-* 


1.701 X 10- 1 


0.3 


3.05 X 10- 1 


5.94 X lO" 1 


1.09 X 10" 1 


2.84 X 10"' 


0.325 


3.95 X lO" 1 


8.4 X 10-' 


1.57 X 10" 1 


4.5 X 10-' 


0.35 


5.0 X 10" 1 


1.13 


2.1 X lO" 1 


6.5 X 10-> 
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Table X — Total and scattering cross sections at 18.1 GHz 
with a = 90° for different drop sizes 



d(cm) 


Qi(cm)* 


Qp(cm)' 


Ql(cm)« 


Q?(cm)« 


0.025 


3.5118 X 10" 6 


3.6944 X 10" 5 


3.7881 X 10"' 


4.0173 X 10-' 


0.05 


5.7214 X 10-* 


6.2236 X 10- 4 


2.4119 X 10-* 


2.7242 X 10- 5 


0.075 


4.4314 X 10- 3 


4.9810 X 10~ 3 


2.8514 X 10-« 


3.4542 X 10"« 


0.1 


1.9770 X 10" 2 


2.2471 X 10- 2 


1.7668 X 10-» 


2.3060 X 10" 3 


0.125 


4.4052 X 10" 2 


5.0478 X 10- 2 


7.1078 X 10" 3 


1.0013 X 10-» 


0.15 


7.4684 X 10" 2 


9.4745 X 10- 2 


1.9966 X 10-» 


3.0950 X 10" 2 


0.175 


1.2199 X 10-' 


1.7421 X 10-' 


4.3974 X 10" 2 


7.4834 X 10 -2 


0.2 


1.8976 X 10" 1 


2.9648 X 10-' 


8.312 X 10" 2 


1.5229 X 10-» 


0.225 


2.782 X 10-' 


4.572 X 10"' 


1.377 X 10-' 


2.640 X 10- 1 


0.25 


3.766 X 10-' 


6.273 X 10"' 


1.999 X 10"' 


3.888 X 10-' 


0.275 


4.714 X 10-> 


7.782 X 10-' 


2.598 X 10" 1 


5.028 X 10- 1 


0.3 


5.56 X 10-» 


9.08 X 10- 1 


3.12 X 10"» 


6.00 X 10-' 


0.325 


6.32 X 10-' 


1.031 


3.54 X 10~ l 


6.89 X 10- 1 


0.35 


7.0 X 10"' 


1.17 


3.9 X 10"' 


7.8 X 10"' 



olated least-squares results by at most \ percent, which is reasonably 
consistent with the error to be expected from the extrapolation. 
Oguchi 34 has since redone his calculations, and he agrees with our 
perturbation results. The same good agreement was obtained between 
the extrapolated and perturbation values of £F(0), S[\>(0), Ql w , and 
Ql Ill) . 

After the least-squares-fitting program had been checked in the 
above manner, we carried out calculations for oblate spheroidal rain- 
drops corresponding to (63), with a and v given by (64). PI ere a is the 



Table XI — Total and scattering cross sections at 30 GHz 
with a = 90° for different drop sizes 



a(cm) 


Q'(cm) 2 


Q« n (om)« 


Qj(cm)* 


0?(cm)« 


0.025 


1.0986 X 10-< 


1.1534 X 10-» 


2.8823 X 10-" 


3.0581 X 10-" 


0.05 


2.2241 X 10"' 


2.4147 X 10" s 


1.9542 X 10- 4 


2.2197 X 10-* 


0.075 


1.4573 X 10- 2 


1.6256 X 10~ 2 


2.5201 X 10- 3 


3.1009 X 10" 3 


0.1 


4.2701 X 10- 2 


5.1454 X 10- 2 


1.4263 X 10" 2 


1.9185 X 10" 2 


0.125 


9.4713 X 10- 2 


1.2407 X 10-' 


4.3818 X 10- 2 


6.2849 X 10- 2 


0.15 


1.6465 X 10-' 


2.1977 X 10-' 


8.9103 X 10" 2 


1.2903 X 10-' 


0.175 


2.3469 X 10"' 


3.0728 X 10-' 


1.3595 X 10-' 


1.9338 X 10"' 


0.2 


2.9521 X 10-' 


3.8199 X lO-i 


1.7408 X 10" 1 


2.4628 X 10"' 


0.225 


3.5403 X 10"' 


4.6666 X 10-' 


2.0731 X 10-' 


3.0183 X 10-' 


0.25 


4.236 X 10"' 


5.800 X 10-' 


2.450 X 10-' 


3.766 X 10-' 


0.275 


5.095 X 10-» 


7.146 X 10" 1 


2.938 X 10-' 


4.710 X 10-' 


0.3 


6.05 X 10-> 


8.47 X 10-' 


3.51 X 10-' 


5.69 X lO-i 


0.325 


7.01 X 10-' 


9.74 X 10" 1 


4.08 X lO-i 


6.62 X lO-i 


0.35 


7.9 X 10- 1 


1.11 


4.6 X 10-' 


7.6 X lO-i 
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Table XII — Total and scattering cross sections at 30 GHz 
with a = 70° for different drop sizes 



a (cm) 


Qi(cm)» 


QP(cm)» 


Qj(cm) 2 


Q, n (cm)* 


0.025 


1.1039 X 10-* 


1.1523 X 10-* 


2.9029 X 10" 6 


3.0582 X 10- 6 


0.05 


2.2367 X 10" 3 


2.4049 X 10- 3 


1.9851 X 10" 4 


2.2196 X 10" 4 


0.075 


1.4701 X 10" 2 


1.6188 X 10- 2 


2.5867 X 10- 3 


3.1000 X 10- 3 


0.1 


4.3613 X 10~ 2 


5.1354 X 10" 2 


1.4859 X 10- 2 


1.9213 X lO" 2 


0.125 


9.8252 X 10-* 


1.2427 X 10-' 


4.6292 X lO" 2 


6.3156 X 10" 2 


0.15 


1.7202 X 10-' 


2.2095 X 10" 1 


9.4694 X 10" 2 


1.3017 X 10- 1 


0.175 


2.4552 X 10" 1 


3.1006 X 10-' 


1.4482 X lO" 1 


1.9599 X lO" 1 


0.2 


3.0950 X 10- 1 


3.8653 X 10- 1 


1.8634 X 10"' 


2.5071 X 10" 1 


0.225 


3.7305 X 10-' 


4.7231 X 10" 1 


2.2389 X 10" 1 


3.0787 X 10" 1 


0.25 


4.496 X lO" 1 


5.861 X 10- 1 


2.676 X 10" 1 


3.837 X lO" 1 


0.275 


5.442 X 10-' 


7.223 X 10-' 


3.240 X 10"' 


4.796 X lO" 1 


0.3 


6.50 X 10-' 


8.59 X 10-' 


3.90 X 10" 1 


5.81 X lO" 1 


0.325 


7.57 X 10" 1 


9.91 X 10" 1 


4.58 X 10" 1 


6.79 X 10- 1 


0.35 


8.6 X 10-' 


1.13 


5.2 X 10" 1 


7.8 X lO" 1 



radius (in centimeters) of the equivolumic spherical drop, and the cal- 
culations were done for d = 0.025(0.025)0.35. The corresponding values 
of a and v are given in Table I. The values taken for the wavelength 
\ = 2ir/k were (in centimeters) 7.5, 2.727, 1.6575, and 1.0, correspond- 
ing approximately to frequencies of 4, 11, 18.1, and 30 GHz. At 20 °C, 
the refractive indices N = 7.884 + 2.184z at 11 GHz, N = 6.859 
+ 2.716* at 18.1 GHz, and N = 5.581 + 2.848i at 30 GHz were ob- 
tained from an elaborate fitting equation in a recently published survey 8 
of available measured data. Since the calculations at 4 GHz were made 
at an earlier date, the value N = 8.77 -f- 0.915i, taken from the older 

Table XIII — Total and scattering cross sections at 30 GHz 
with a = 50° for different drop sizes 



a (cm) 


Q«(cm) 2 


QP(cm) 2 


Qj(cm)* 


Q^cm)' 


0.025 


1.1173 X lO" 4 


1.1495 X 10- 4 


2.9552 X lO" 6 


3.0584 X 10- 6 


0.05 


2.2684 X lO' 3 


2.3803 X 10" 3 


2.0634 X lO" 4 


2.2193 X 10-« 


0.075 


1.5025 X lO" 2 


1.6015 X 10" 2 


2.7560 X 10" 3 


3.0978 X 10" 3 


0.1 


4.5936 X lO" 2 


5.1102 X lO" 2 


1.6378 X 10- 2 


1.9286 X lO" 2 


0.125 


1.0732 X lO" 1 


1.2476 X 10-' 


5.2632 X lO" 2 


6.3941 X 10- 2 


0.15 


1.9102 X 10-' 


2.2399 X 10"' 


1.0914 X 10" 1 


1.3308 X lO" 1 


0.175 
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Fig. 3 — Comparison of perturbation approximations to the least-squares-fitting 
value of Re <Si(0) at 4 GHz with a = 90° as a function of drop size. 

literature, was used, rather than N = 8.78 + 0.977*. The angle of inci- 
dence a was taken to be 90° at 4, 11, and 18.1 GHz, while at 30 GHz 
the calculations were done for a = 70° and a = 50° also. 

The calculated values of the forward scattering amplitudes S\(0) 
and Sii (0) are given in Tables II to VII, and those of the total cross 
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sections Q\ and Q} 1 and the scattering cross sections Q] and QV are given 
in Tables VIII to XIII, all rounded in the last significant figure. The 
values of the absorption cross sections Ql and Ql l follow from (38). 
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Fig. 4 — Comparison of perturbation approximations to the least-squares-fitting 
value of Re S\\ (0) at 4 GHz with a = 90° as a function of drop size. 
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The accuracy of the least-squares fit of the boundary conditions de- 
creases with increasing drop size, because of the increase in eccentricity, 
which is why fewer significant figures are given in the tables for the 
larger drop sizes. Except for the smaller drop sizes, for which the results 
could be given more accurately, the number of significant figures 
reflects the degree of convergence of the results, as evidenced by in- 
creasing the upper limit of n in the least-squares fit by 2 and by 4. The 
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Fig. 5 — Comparison of perturbation approximations to the least-squares-fitting 
lues of Re Si (0) and Re S u (0) at 11 GHz with a = 90° as a function of drop size. 
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accuracy of the far-field results is generally at least one order of mag- 
nitude greater than that of the fit of the boundary conditions for the 
reasons discussed in Section V. We note that the largest drops occur 
only at the heaviest rain rates, 7 and then only a small percentage of 
them, so that the lower accuracy of the results for these drops is ac- 
ceptable when summing over the drop size distribution. 

The number of terms required to obtain the desired accuracy for the 
far-field quantities and to adequately satisfy the boundary conditions 
increases with both drop size and with frequency. At 4 GHz, it was 
found that max m = 4 and max n =17 were sufficient for the largest 
drop size. For a = 90° at 30 GHz, it was necessary to take max m = 8 
and max n = 23 for the largest drop size. In this latter case, more than 
half the capacity of the Honeywell 6070 computer was used. In some 
cases, it was found that max n or max m could not be increased without 
causing overflow in some of the subroutines, in particular SBES and 
L2FIT. 
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Fig. 6 — Comparison of perturbation approximations to the least-squares-fitting 
value of Re 5i(0) at 18.1 GHz with a = 90° as a function of drop size. 



988 THE BELL SYSTEM TECHNICAL JOURNAL, JULY-AUGUST 1974 



To check on the advantage of using least-squares fitting (with ap- 
proximately twice as many fitting points as unknown coefficients) 
rather than collocation, we used collocation in several cases at different 
frequencies and for different drop sizes. Our general conclusion is that, 
for the same max m and max n, results may be obtained by collocation 
for the far-field quantities that are almost as accurate as those obtained 
by least-squares fitting. However, there are much larger errors in the 
boundary conditions (in between the fitting points) with collocation 
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Fig. 7 — Comparison of perturbation approximations to the least-squares-fitting 
value of Re Sn(0) at 18.1 GHz with o = 90° as a function of drop size. 



ELECTROMAGNETIC WAVE SCATTERING 989 



than with least-squares fitting. For the larger, more eccentric raindrops, 
some of these errors were of the order of 100 percent, which seem to be 
unacceptable. However, for the smaller raindrops, the errors in the 
boundary conditions with collocation are acceptable. Since the cost 
of carrying out the least-squares fit is less when fewer fitting points are 
used, collocation has the advantage of reducing the cost, although fewer 
terms are required, anyway, to obtain the desired accuracy for the 
smaller raindrops. It is possible that the collocation fit may be improved 
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Fig. 8 — Comparison of perturbation approximations to the least-squares-fitting 
values of Re Si (0) and Re Sn (0) at 30 GHz with a = 90° as a function of drop size. 
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by satisfying the boundary conditions at nonuniformly spaced points, 
but we have not investigated this. 

As a check on the point-matching (collocation) results of Oguchi, 6 
we carried out the least-squares fitting for a = 90°, at 19.3 GHz for 
a = 0.3 and at 34.8 GHz for a = 0.075, 0.15, 0.225, and 0.3, using 
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Fig. 9 — Comparison of perturbation approximations to the least-squares-fitting 
values of Re Si(0) and Re Sh(0) at 30 GHz with a = 50° as a function of drop size. 
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Fig. 10 — Comparison of perturbation approximations to the least-squares-fitting 
values of Im Si(0) and Im Sn(0) at 4 GHz with a = 90° as a function of drop size. 
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Oguchi's relationship 5 



b = { l -u a )> 



(74) 



instead of the first relationship in (1). Our results for the forward scat- 
tering amplitudes are consistent with his point-matching values, but 
they may be given to greater accuracy. Our truncated values for 
a = 0.3 are given below where, in Oguchi's notation, 5 

-5A ., 



f" X 10 3 = 



tSI(O), 



jh x 10 3 = _^^( ), 



(75) 
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Fig. 11 — Comparison of perturbation approximations to the least-squares-fitting 
values of Im -Si (0) and Im Su (0) at 11 GHz with a = 90° as a function of drop size. 
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Fig. 12 — Comparison of perturbation approximations to the least-squares-fitting 
lues of Im Si (0) and Im Sn (0) at 18.1 GHz with a = 90° as a function of drop size. 
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and the vertical lines indicate where Oguchi truncated his results : 



GHz 
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The values taken 34 for the wavelength A (in centimeters) were 1.5533330 
and 0.86135810, corresponding approximately to frequencies of 19.3 
and 34.8 GHz, with refractive indices N = 6.5449188 + 2.8104040i and 
N = 5.0487284 + 2.7948416i, respectively. 
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Fig. 13 — Comparison of perturbation approximations to the least-squares-fitting 
values of Im Si (0) and Im Su (0) at 30 GHz with a = 90° as a function of drop size. 
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The main reason for the greater accuracy of our results is that we 
took larger values of max n than Oguchi, who did the point-matching 
at both frequencies for max n = 12 and 14, with max m = max n. For 
d = 0.3 we took max m = 7 at 19.3 GHz and max m . = 9 at 34.8 GHz, 
which were sufficient, and took max n = 21 at both frequencies for the 
least-squares fitting. We also used collocation for a = 0.3 and max 
n = 12, 14, and 21. For max n = 21 the collocation results differ by 
at most 1 in the last decimal place from the results given above, but 
errors in some boundary conditions were of the order of 10 percent, as 
compared with much less than 1 percent for least-squares fitting. This 
is consistent with our general conclusion discussed earlier in this sec- 
tion. We point out that the raindrops satisfying (74) are less eccentric 
than those satisfying (1), so that the overall errors are correspondingly 
smaller. For collocation with max n = 14, some errors in the boundary 
conditions were close to 100 percent, which explains why, with point- 
matching, Oguchi did not give any significant figures for Re f h at 
34.8 GHz, for either a = 0.3 or a = 0.325. 

Although Oguchi 5 gives four significant figures for f v and f h at 34.8 
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Fig. 14 — Comparison of perturbation approximations to the least-squares-fitting 
values of Im Si (0) and Im Sn (0) at 30 GHz with a = 50° as a function of drop size. 



996 THE BELL SYSTEM TECHNICAL JOURNAL, JULY-AUGUST 1974 



GHz for all drop sizes corresponding to his solution in terms of sphe- 
roidal wave functions (with modal sums truncated at 9), these values 
are not consistent with his point-matching ones for the larger drop 
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Fig. 15 — Comparison of perturbation approximations to the least-squares-fitting 
value of Q} 1 — Q\ at 4 GHz with a = 90° as a function of drop size. 
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sizes. For a = 0.3, he gives f h X 10 3 = 0.06470 - 4.709i, which is in 
fact quite close to our value, but he also gives /" X 10 3 = 0.9100 
— 3.726i, with real part differing by almost 1 percent from our value. 
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Fig. 16 — Comparison of perturbation approximations to the least-squares-fitting 
value of Q} 1 — Q, at 11 GHz with a = 90° as a function of drop size. 
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VIII. PERTURBATION RESULTS 

In this final section, we compare three sets of first-order perturbation 
results with those obtained by least-squares fitting. The comparisons 
are made graphically in Figs. 3 to 23, since this is much more revealing 
than tabulating the results. The solid curves correspond to least- 
squares fitting and the dashed curves to perturbation about the in- 
scribed sphere of radius o, corresponding to the expansion in (63), with 
perturbation parameter v given by (64). The circles and dots corre- 
spond to perturbation about the equivolumic sphere of radius a, with 
perturbation parameters v = '2a and v, respectively, corresponding to 
the expansions in (65). The dots have been omitted in those cases in 
which they would lie very close to the corresponding circle or solid 
curve. Comparisons are made for a = 90° at 4, 11, 18.1, and 30 GHz 
and for a = 50° at 30 GHz. 

The real parts of the forward scattering amplitudes Si (0) and Sn (0) 
and the first-order approximations to these quantities are depicted in 
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Fig. 17 — Comparison of perturbation approxi 
,lue of Q 11 - Q] at 18.1 GHz with a = 90° as : 
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Fig. 18 — Comparison of perturbation approximations to the least-squares-fitting 
value of Q l \ — Q\ at 30 GHz with a = 90° and a = 50° as a function of drop size. 

Figs. 3 to 9, while the imaginary parts are depicted in Figs. 10 to 14. 
It should be noted that a logarithmic scale has been used in Figs. 3 
and 4. Thus, at 4 GHz the first-order approximation to the real part of 
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Fig. 19 — Comparison of perturbation approximations to the least-squares-fitting 
value of Im[Si(0) - <S'n(0)J at 4 GHz with a = 90° as a function of drop size. 

Sn(O), obtained by perturbing about the inscribed sphere, is in error 
by an order of magnitude for the largest drop size. It is seen that the 
best overall approximation is obtained by perturbing about the 
equivolumic sphere with perturbation parameter i> = 2d, and in most 
cases there is a significant improvement over the approximation ob- 
tained by perturbing about the inscribed sphere as Oguchi 4 did. The 
second best overall approximation is obtained by perturbing about the 
equivolumic sphere with perturbation parameter v = o(2 — a), and is 
generally much better than the approximation obtained by perturbing 
about the inscribed sphere. The above ordering of the three sets of 
perturbation results is consistent with the order of the geometrical 
errors in the corresponding approximations in (63) and (65) to the 
oblate spheroid. 

Although the comparison is not depicted for some of the smallest 
drop sizes, all three approximations are good for these, since the ec- 
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Fig. 20 — Comparison of perturbation approximations to the least-squares-fitting 
value of Im[*Si(0) — Sn(0)J at 11 GHz with a = 90° as a function of drop size. 

centricity is small. On the other hand, the approximations obtained 
by perturbing about the equivolumic sphere are remarkably good for 
the largest drop sizes, in view of the fact that neither the eccentricity 
nor the perturbation parameter is small. In particular, these approxi- 
mations to the imaginary part of $n(0), depicted in Figs. 10 to 14, are 
quite impressive. It is not too surprising that perturbing about the 
inscribed sphere leads to poor results for the larger drop sizes in the 
second polarization. The first-order approximations to the scattering 
cross sections Ql and Ql 1 are very similar to those depicted in Figs. 3 
to 9 for the real parts of Si (0) and Su (0), which are related to the total 
cross sections Q] and Q} 1 by (39). 
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Fig. 21 — Comparison of perturbation approximations to the least-squares-fitting 
value of Im[Si(0) - .Sn(0)] at 18.1 GHz with a = 90° as a function of drop size. 

For purposes of comparison, the values of the forward scattering 
amplitude S(0) and the total and scattering cross sections Q, and Q s 
for the equivolumic spherical drops are given in Tables XIV to XVII. 
These quantities do not depend on the polarization of the incident 
wave or on the angle of incidence a. As is seen from Tables VIII to XIII, 
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Fig. 22 — Comparison of perturbation approximations to the least-squares-fitting 
value of Im[<Si(0) - Sn(0)] at 30 GHz with a = 90° as a function of drop size. 
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Fig. 23 — Comparison of perturbation approximations to the least-squares-fitting 
value of ImCSi^O) - 5n(0)J at 30 GHz with a = 50° as a function of drop size. 

the value of Q t lies between the corresponding values for the two polar- 
izations for the oblate spheroidal drop of the same size and similarly 
for the value of Q s . Although this happens to be true at 30 GHz for 
a = 90°, 70°, and 50°, these relations should not be expected to hold 
for all values of a, since for a = 0° the cross sections are independent 
of the polarization because of the axial symmetry of the oblate sphe- 
roidal drop. We have verified that Q t ^ Q\ = Q? and Q t ^ Q] = QV 
for a = 0° at 11 GHz, for a = 025, 0.1, and 0.175. 

The rain-induced differential attenuation and differential phase 
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Table XIV — Forward scattering amplitude and total and 

scattering cross sections for the equivolumic spherical 

drop at 4 GHz for different drop sizes 
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shift are obtained 1 by summing the real and imaginary parts of $11 (0) 
— Si (0) over the Laws and Parsons drop-size distribution. 7 In a recent 
short note, 9 the three first-order perturbation approximations have been 
compared to the least-squares fitting results for the differential atten- 
uation and differential phase shift at several different rain rates. The 
same ordering of the overall closeness of the three approximations 
holds for these quantities. Since the perturbation results are obtained 
quite inexpensively whereas the least-squares-fitting procedure is 
very costly, approximations to the differential attenuation and differ- 
ential phase shift at frequencies up to 100 GHz were obtained by per- 

Table XV — Forward scattering amplitude and total and 

scattering cross sections for the equivolumic spherical 

drop at 11 GHz for different drop sizes 



a (cm) 


5(0) 


Q<(cm) 2 


Q.(cm) 2 


0.025 


4.9868 X 10-* 


— 1.8550 X 10-4 


1.1804 X 10" 6 


5.3599 X lO"" 8 


0.05 


6.2648 X 10 -6 


- 1.5238 X 10- 3 ? 


1.4829 X lO" 4 


3.4733 X lO" 6 


0.075 


3.8512 X lO" 4 - 


- 5.3789 X lO" 3 * 


9.1163 X 10" 4 


4.0552 X 10-* 


0.1 


1.7992 X 10" 3 - 


- 1.3547 X 10- 2 i 


4.2588 X 10- 3 


2.3883 X 10~ 4 


0.125 


7.1756 X lO" 3 - 


- 2.7813 X lO" 2 ? 


1.6985 X 10" 2 


1.0004 X lO" 3 


0.15 


2.2023 X lO" 2 - 


- 4.4760 X 10- 2 j 


5.2132 X 10" 2 


3.4809 X 10" 8 


0.175 


3.9051 X 10~ 2 - 


- 5.5463 X 10" 2 * 


9.2438 X 10~ 2 


9.7067 X lO' 3 


0.2 


5.0153 X 10- 2 - 


— 7.3383 X 10- 2 i 


1.1872 X lO" 1 


2.1692 X lO" 2 


0.225 


6.8317 X lO" 2 - 


- 1.0343 X 10-4 


1.6171 X 10- 1 


4.2294 X 10" 2 


0.25 


9.6583 X lO" 2 - 


— 1.3932 X lO -1 ? 


2.2862 X lO' 1 


7.4229 X lO" 2 


0.275 


1.3508 X lO" 1 - 


- 1.8178 X lO" 1 ? 


3.1976 X lO" 1 


1.2397 X lO" 1 


0.3 


1.8972 X 10-' - 


- 2.2945 X 10-4 


4.4910 X lO" 1 


2.0165 X lO" 1 


0.325 


2.6372 X lO" 1 - 


- 2.7484 X 10-4 


6.2427 X 10" 1 


3.1470 X 10" 1 


0.35 


3.5439 X 10" 1 - 


-3.1001 X 10-4 


8.3887 X 10" 1 


4.6135 X 10" 1 
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Table XVI — Forward scattering amplitude and total and 

scattering cross sections for the equivolumic spherical 

drop at 18.1 GHz for different drop sizes 



d(cm) 


5(0) 


Q ( (cm) 2 


Q,(cm) 2 


0.025 


4.1444 X 10- 5 


— 8.3170 X 10- 4 ?: 


3.6243 X 10- 5 


3.9388 X lO" 7 


0.05 


6.8431 X 10" 4 


— 6.9938 X lO' 3 / 


5.9843 X 10- 4 


2.6138 X lO" 5 


0.075 


5.3548 X 10- 3 


— 2.4905 X lO" 2 / 


4.6827 X lO -3 


3.2305 X 10" 4 


0.1 


2.4004 X 10" 2 


— 5.561 1 X 10-*t 


2.0991 X 10~ 2 


2.0935 X lO" 3 


0.125 


5.4019 X 10~ 2 


— 9.1780 X 10-*i 


4.7240 X lO" 2 


8.8579 X lO" 3 


0.15 


9.7706 X 10- 2 


— 1.4719 X 10-'/ 


8.5444 X 10" 2 


2.6617 X 10- 2 


0.175 


1.7343 X 10" 1 


— 2.1684 X lO -1 / 


1.5166 X 10-' 


6.2646 X 10~ 2 


0.2 


2.8868 X 10- 1 


— 2.8757 X 10->i 


2.5245 X 10" 1 


1.2552 X 10" 1 


0.225 


4.4398 X 10-' 


— 3.3858 X lO-'i 


3.8826 X 10"' 


2.1846 X 10-' 


0.25 


6.1859 X 10-' 


— 3.5216 X 10-»i 


5.4095 X 10"' 


3.2834 X 10-' 


0.275 


7.8546 X 10-' 


— 3.3543 X lO-'t 


6.8688 X 10-' 


4.3596 X lO" 1 


0.3 


9.3560 X 10-' 


— 3.0993 X i0~H 


8.1818 X lO"' 


5.3377 X lO' 1 


0.325 


1.0751 


— 2.9126 X lO -1 / 


9.4019 X 10" 1 


6.2372 X lO" 1 


0.35 


1.2136 


— 2.8673 X lO-'t 


1.0613 


7.1027 X lO" 1 



turbing about the equivolumic sphere. However, the results may be 
less reliable at the higher frequencies, particularly at the heavier rain 
rates. 18 

The difference Q" — Q\, which is related to the real part of Sn(0) 
— Si(0) by (39), is depicted in Figs. 15 to 18, and the imaginary part 
of Si(0) - flu (0) is depicted in Figs. 19 to 23. We note that, although 
extra first-order correction terms arise in the expansions about the 
equivolumic sphere given in (65), they correspond to a constant change 
in the radius of the drop. Hence, the corresponding increments in the 
forward scattering amplitudes are the same for both polarizations, and 

Table XVII — Forward scattering amplitude and total and 

scattering cross sections for the equivolumic spherical 

drop at 30 GHz for different drop sizes 



a (cm) 


5(0) 


Q«(cm) 2 


Q.(cm) 2 


0.025 


3.5549 X 10- 


4 — 3.8212 X 10- s i 


1.1316 X 10 ' 


2.9980 X 10- fi 


0.05 


7.2950 X 10- 


1 — 3.2466 X 10-H* 


2.3221 X lO" 3 


2.1254 X lO" 4 


0.075 


4.8577 X 10" 


2 — 1.0204 X lO-'i' 


1.5463 X 10- 2 


2.8859 X 10-' 


0.1 


1.5024 X 10- 


1 — 2.0855 X 10-' i 


4.7823 X 10- 2 


1.7391 X 10- 2 


0.125 


3.5560 X 10- 


' — 3.1825 X 10-'/ 


1.1319 X 10-' 


5.6317 X lO" 2 


0.15 


6.3545 X 10" 


• — 3.5742 X 10-V 


2.0227 X 10-' 


1.1723 X 10-' 


0.175 


9.0831 X 10- 


' — 3.2571 X 10-'/ 


2.8913 X 10-' 


1.8042 X 10-' 


0.2 


1.1444 


— 2.9371 X 10"'/ 


3.6426 X 10-' 


2.3403 X 10"' 


0.225 


1.3868 


— 3.0778 X 10-'/ 


4.4142 X 10-' 


2.8578 X 10-' 


0.25 


1.6859 


— 3.5609 X 10-' 1 


5.3662 X 10-' 


3.4813 X 10-' 


0.275 


2.0559 


— 3.9425 X 10-' i 


6.5441 X lO' 1 


4.2750 X 10-' 


0.3 


2.4657 


— 3.9064 X 10-'/' 


7.8486 X 10-' 


5.1929 X 10"' 


0.325 


2.8763 


— 3.5872 X 10"'/ 


9.1557 X 10-' 


6.1372 X 10-' 


0.35 


3.2831 


— 3.3474 X lO" 1 / 


1.0450 


7.0680 X 10"' 
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therefore do not affect the difference £n(0) - Si(0). Figures 15 to 23 
show that the approximations to the differential quantities obtained 
by perturbing about the equivolumic sphere with perturbation param- 
eter v = 2d are overall remarkably close to the least-squares-fitting 
results and far better than the approximations obtained by perturbing 
about the inscribed sphere. 
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APPENDIX A 

We first derive the expansion of the incident plane wave in a Fourier 
series in the azimuthal angle <p, as given by (20). The unit vectors in 
Cartesian coordinates are given in terms of those in spherical coordi- 
nates by 

i = sin cos <pii + cos cos <pi 2 — sin <pi 3 , 

j = sin sin <pii + cos sin <pi 2 + cos <pia, 

and 

k = cos 0ii — sin 0i 2 . C 76 ) 

Also, we have 

x sin a + z cos a = r(sin a sin cos <p + cos a cos 0). (77) 

But 36 for integer values of p, 

— / e-'p^exp^cos (p)d<p = i p J P (t), (78) 

2,-K Jo 

where J v denotes the regular Bessel function (of the first kind) of 
order p. 
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It follows, from (5), (6), (20) to (22), and (76) to (78), using the 
recurrence relations for the Bessel functions, 37 that 

fm(r, 6) = i m exp (ik r cos a cos 0) J m (k r sin a sin 0) sin a 

X (sin 0i 2 — cos 0ii) — iJ' m {kor sin a sin 0) cos a (sin 0ii + cos 0i 2 ) 

■ mJ m (k r sin a sin 0) . "1 ,_ . 

H r^ — = =— a — cos ai 3 (79) 

k r sm a sin J 

and 

/ fl v / - 7 M T m J™ ( k <> r sm a sm 0) 
gm(r, 0) = —i m exp (ik r cos a cos 0) ; = ; — - — 

b L V sin a sin 9 

X (sin 0ii + cos 0i 2 ) + il' m (kor sin a sin 0)i 3 , (80) 



where, as before, the prime denotes derivative with respect to the 
argument. 

From (8) to (11), (13), (14), and (20), the boundary conditions (16) 
and (17), when multiplied by e~ imv and integrated with respect to <p 
from to 2ir, lead to the equations 

e m3 (R, 0) + E flmn/i» («oft) jt 

rv^|m| Off 

715^0 

- E ^^^ + Ai 1) '(*^)l-^ S Pl" l (oo8fl) 

n§ | m | L koti J sin0 



1^0 

= E C mn j n {kyR) 
715*0 



dPj, m| (COS 0) 

d0 



n&|m| L *1** J 

715*0 



vm 
sin 



Pl m| (cos0) (81) 



and, using (4), 

^h nZ {R,6)+ E h mn h£\k»R) 

#0 ng|m| 

715*0 



dPJ, m| (COS 0) 
^0 



- £ a„I*^ + A4»'(W)l--^P!,»l(oo S 9) 

7iS|m| L Kott J Sln " 

dPJ, m| (cos 0) 



715*0 



= iV Z d mn j n {krR) 

71 ^ |77l | 
715*0 

-N E c mn 



n 2.\m\ 

715*0 



^» + ii < W8 ,].^rt-(„fl. (82) 

We are using the notations e mj - = e m i, and h m j = h m -iy. 
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Similarly, the boundary conditions (18) and (19) lead to the 
equations 

e mi (R, B)+± d § e ml (R, 6) - £ a mn h™ (koR) ~ Pi"' (cos 6) 
K ao n > i m i sin u 



n2|m| 



(k R) 



\\ h™(k l 
in \[ k R 



+ K iy (k R) 



n § |m| 



IdPT 



(cos 0) 



d0 



P d0 fc P 



= - E c mn j„(* 1 fl)--^PW(cos0) 

nS|m| Sni 



(cos 0) 



d0 



71^0 



+ n(n + l)gi,(^fi) pi „ (co8g) 



P d0 A^P 



(83) 



and 



fc>>Mo 

A; 



[*-(«,«+ g^ *-i(B,«)] 



Z??7 



(cos 0) 



d0 



- L 6«^»(AoB)-^Prl(cosff) 

n & | m | Sin 

nS|m| l|_ fc 0« J 

n?50 

nS|m| SU10 

dPi"* 1 (COS 0) 



np*0 



dd 



n(n + 1) <2P jnjkjR) plm{ f fl v 

+ r ~di~kjr^ n {cos0) 



(84) 



APPENDIX B 

We consider here the case in which the raindrop is symmetrical 
about the plane = ir/2 so that 

R( T - 6) = R(d), OgOg ir/2. (85) 

Let a = (ir - a). Then, from (21), (79), and (80), it follows that, cor- 
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responding to <5, 

&»(#,*) = -<Z*(R,t-0), (86) 

hUR,o) = h l m3 (R,*-e), (87) 

[& 2 (fl,0)+]^^i(fl,0)] 

= [eL(/2, tt - d) + Yd{*-d) Al{R ' * ~ $) ]> (88) 
and 

= - [h l m2 (R, w - 0) + ^ ^£^ ttito »-*)]■ (89) 

But" 

Pj, w| (-cos0) = (-l)» + i"*'Pi M| (cos0). (90) 

It follows from (81) to (90) that 

«&. + (-l)» + '"' +1 C - 0, cL+ (-l) n+|m|+, ^ = (91) 
and 

ft, + (-D n+|m| bL = 0, A. + (-l) n+|m| ^n = 0. (92) 

For a = 7T/2, we have & = a, and hence we obtain the relationships 
in (25). 
Similarly, from (22), (79), and (80), 

t&(R,0) = e&iR, r - 6), (93) 

h'MR, 0) = -h%s(R, w-0), (94) 

[fi(B,»)+5§fi(«,«] 

= - [<&(R, * - 6) + i ^T^ AW *-*)]» (95) 
and 

= [*&(*, * - 0) + i ^J^ *&<*, »-«)]• (96) 
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It follows, from (81) to (85), (90), and (93) to (96), that 

<&, + (-1)»+'*I<J£. = 0, (& + (-1)*H1£, = (97) 

and 

b l L + (-l)» + "»'+ l $Si, = 0, dSn + (-1)«*»-U* = 0. (98) 

For a = tt/2 we obtain the relationships in (26). 

For a 7* 7r/2 we may consider the sum and the difference of the 
boundary conditions corresponding to a and to & = (ir — a). Then 
the sums (a mn + d mn ), (b mn + b mn ), (c mn + t nn ), and (d mn + d nn ) 
and the differences (a mn — d mn ), (b mn — b mn ), (c mn — £ mn ), and 
(d m „ — d mn ) may be determined separately, and these sums and 
differences vanish for alternate values of n, depending on the pol- 
arization. 

APPENDIX C 

We first consider the calculation of the scattered energy W„ which 
is defined by (35), by letting r — ><x . From (29) it follows that 

W. = lim \^- C [' (|#J|»+ \Ei\*) sine dd d<p\ ■ (99) 

r _oo [ Z0)fi Jo Jo | 

But from (28), 

oikor <o r rn . . 

E '*~-jrr £ s (-0" U-.-^^'Ccostf) 

/C r m =-«, nfc|m| L Sm ^ 

+ t „^o^]^ (100 ) 



and 



pikar =o I 

rt-0' m= — co n&|m| L 



dPi ro| (cos e) 



dd 

n*0 



+ 6mn • -^-, Pi" 1 (cos 0) |e<"MP. (101) 

Sill (7 



I le""*'. 



Substituting (100) and (101) into (99), the integration with respect 
to <p is straightforward. The integration with respect to 6 readily follows 
with the help of the identities 



ml |Pi"l(co 8 ») (U>i " l(c08 ' ) 



/;[pj-i(, 



dd 

, dPl" 

dd 
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+ < "T£ e "* ) f.-'(eo.»)}»-0 (102) 



3d 



and 



/;[ 



dPl^pst) dP!r> (cosy w? pJ „, ,„, os i in m 

dd dd sin 2 J 

2n{n + \ ){n + \m\)\ Mn „v 

= (2n + l)(n- |m|)! 5 "" (1 ° 3) 



where 5i„ denotes the Kronecker delta, i.e., Si„ = 1 for J = n, and 
otherwise. Thus, the expression for W, given in (36) is obtained. 

We remark that there is no need to let r — » oo to obtain this expression 
for W,. The same result follows from (35) by using the expressions (10) 
and (11) for the scattered field, wherein MJS(fc ) and N£l(k ) are 
defined by (8) and (9), with z„(k r) = A£°(fc<r). The dependence on r 
is found to vanish, as is to be expected, in view of the Wronskian rela- 
tionship 39 

jn{k<r)yL{k(/r) - yn{k r)jn{k r) = 7^y 2 " ( 104 ) 

We also remark that the expression in (36) holds quite generally, e.g., 
for scattering from nonaxisymmetric raindrops, since at this point we 
have made no use of the properties of the coefficients a mn and b mn . 

We next consider the calculation of the total energy W t , which is 
defined by (41). We begin by allowing for a general incident field, 
given by 

E' = - t E [4«Mffi(*«) + B.JttJK*.)] (105) 

m = — «> n& |m| 

and 

H< - *i £ £ CA*nNS>(*o) + B^MffiCfci)! (106) 

W/Io m = - « ng|m| 

where the superscript 1 indicates that z n (k r) = j n (k r) in (8) and (9). 
The calculation of W t is similar to that of W. and it is found, after 
some reductions, that 

_ -2tt p - ^ n(n + l)(n + |m|)l 
Wt o,Mofco .±l-.^r«| (2n + l)(n- |m|)l 

X(a m Xn + t ra X). (107) 

We now consider the incident electric field given by (42), where 
Ei and Ek are given by (5) and (6), respectively. Then (43) holds and, 
from (10) and (11), 

a mn = (&„ + oS., 6 mn = bL + bg» (108) 
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From (5) and (6) and the expansions in (116) and (117), it follows 
from (105) to (107) that 

W t = -^Re f E (-i)"- 1 

nyiQ 

( r^r w „■_, , N , , dP l n m * (cos a) 1 

x K*--a= rt (cosa) + 6 - — * J 

+ ig;,[a. / P """. (c0Sa) + 6... J2-rt-l(co.«)l} ■ (109) 
L cia sin a J ] 

The relations (39) now follow from (33), (34), (40), and (108) by 
setting first E u = and second E\ = 0, in (109). We also note that, 
from (5), (6), (28), (30), (42), and (109), 






.. 1 -;fcore- i * or (E , ')*-E a U-«, . 
hm ' 



exp[— ik (x sin a + z cos 



§i)]- < 110 > 



We remark that both (109) and (110) hold, subject to (42), (43), and 
(108), for scattering from generally shaped raindrops. 

Finally, we consider the particular case of an axisymmetric raindrop 
given by r = R{d), so that (23) and (24) hold. Thus, 

{a} mn aZ + b&Jfi) - - {a'-mnO^n + hLnjKn). (Ill) 

Hence, from (36) and (108), 

W. = W\ + Wf. (112) 

Then, from (37), with 

Qs = yr, f ( H3) 

ko(EiEj + EuEii) 

we obtain (44). Also, from (23), (24), (108), and (109), 

W t = W\ + W?. (114) 

Hence, from (40), with 

= 2tuAt TFt M15) 

ko(EiEl + EnEii) 

we obtain (45). 

APPENDIX D 

We first give the expansions for the incident wave in terms of spheri- 
cal vector wave functions. 419 It is found that 
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(cos ai — sin ak) exp [ik (x sin a + z cos a)] 

- - i n+l {2?i + l)(n- Im))! 

OT A M „= 2 rm| n(n+ l)(n+ |m|)l 



X [ -^- Pi- 1 (cos «)M<» (Ao) + dP " W(cosa) Nfli (Ao) 1 (H6) 
L sin a da J 

and 

j exp\jik (x sin a + z cos a)] 

- - i»(2n + l)(n - |m|)I 

m il.nir»| n(n+ 1) (» + |m|)! 

X [ dPj,OT ' (C0Sa) Mffl (ifco) + -£- Pi" 1 (cos a)Nfli(A.)l ■ (117) 
[ da sin a J 

Expressions for the quantities e m (r, 0) and h m (r,d), defined in (20), 
then follow from (5), (6), (8), and (9). Thus, we may now consider the 
boundary conditions (81) to (84). 

We first multiply (81) by cZPJ m| (cos 6)/dd sin and (83) by im P] m| 
(cos 0) and add, and then multiply (81) by im P} ml (cos B) and (83) by 
dP\ ml (cos 6)/dd sin 6 and subtract, and integrate both these equations 
with respect to 6 from to it. In the zero-order approximation corre- 
sponding to v = in (47), this leads, with the help of (102) and (103), 
to simultaneous linear equations for a$ and c$$. Similarly, multiplying 
(82) by dPj ro| (cos0)/d0sin0 and (84) by im P\ ml (cos d) and adding, 
and multiplying (82) by im P] m| (cos0) and (84) by dP] m| (cos d)/dd 
sin 6 and subtracting, and integrating both these equations with re- 
spect to 6 from to ir, we obtain simultaneous linear equations for the 
zero-order coefficients b§$ and d$. The solution of these two pairs of 
simultaneous equations leads to the relations (50) and (51), where the 
quantities a mn and /3 mn depend on the polarization, as given by (52) 
and (53). It remains to give the expressions for the quantities a„, 6„, c„, 
and d n occurring in (50) and (51). 

We define 

F.(€) = [ ^ + A. m '<0], Gn(t) = [ ^ + *(«)] • (H8) 

Then, with p = koa, it is found that 

lJ n (p)NG n (N P ) - j n (Np)G n ( P n (nq) 

ih^(p)NG n (N P ) - jn(Np)F n (p)l ' K ' 

h Un(p)Gn(Np)-Njn(Np)G»(f)l ^x 

" ~ W(p)G n (Np) - Nj n (Np)F n (p)l ' K 



a„ = 
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and 

en = (i/p 2 )K ) ( P )NG n (Np) - j n (Np)F n (p)2- 1 , (121) 

d„ = (i/p i )W! ) (p)G n (Np) - Nj n (Np)F n ( P n- 1 . (122) 

In obtaining (121) and (122), we have made use of the Wronskian 
relationship 39 

Mp)W(p) ~ h^{ P )jn{p) = (i/p 2 )- (123) 

Next, considering the first-order terms in v in the integrated forms 
of the boundary conditions, and making use of (47) to (49), two pairs 
of simultaneous linear equations are obtained for aj& and c£l and for 
b££ and d£». These equations contain somewhat involved expressions, 
but after considerable reductions they lead to the expressions given in 
(54) to (56), subject to (57) to (62). In particular, use has been made 
of the differential equation satisfied by the spherical Bessel functions, 39 

e*f (e) + 2&U) + ie - i(i + i)XD = o. (124) 

Moreover, from (118) to (124) it follows that 

ji(p) - aih^(p) + ciji(Np) = 0, (125) 

j'M - ath^'ip) + Ncij'i(Np) = 0, (126) 

G\(p) - aiF'i(p) + WaG'tiNp) = (1 - N*)cij t (Np), (127) 

NjM - NbihfHp) + dMNp) = (1 - N*)diji(N P ), (128) 

j'i(p) ~ W(p) + NWij'iiNp) = (N> - DdtGtiNp), (129) 

and 

Np^G'^p) - biF't(p) + NdiG' t (N P )2 

= J(l + l)(l -N*)drfi(N P ). (130) 

We have also used the fact that 

A Pi (cos e) — de — sin e ~dB de - in + 1) *- '* (131) 

where OC^J and d% are given by (60) and (61). This result follows directly 
by integration by parts and use of the differential equation satisfied 
by the associated Legendre functions, 14 

1 d[ . . d PH(cosfl) "1 

s!n-^L sm * de J 

+ \n(n + 1) - £^- e ]pl ml (cos d) = 0. (132) 
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APPENDIX E 

We outline here the calculation of the integrals in (60) to (62) in the 
case 

<n(0) = |sin 2 0. (133) 

It is assumed that I ^ \m\, n ^ \m\, I ^ and n ^ 0. Integration 
by parts of the expression in (62) leads to 

SZ = -m i T P\ m] (cos 0)PJT I (cos 0) sin cos 0d0 
./o 

= -mj 1 xP\ ml (x)Pl m \(x)dx. (134) 

But," with P|S!-i(a?) = 0, 

(2i+ l)xPj m| (x) = (i - |m| + l)PJTl(z) + (* + |m|)PW(s). (135) 

Substituting (134) into (135), and using the relationship 14 

/_; pj-. W flr< m* = (2w %V-'|li)! { '" (136) 

it follows from (62) that 

Jm -m \ (n - \m\) , (n + \m\ + 1) ^ 1. n«m 

J - = ^TI) [ (2n - 1) 8l - 1 + (2» + 3) 8l - +1 \ (13?) 

Next, from (60) and (133), 

3C5 = g J(« + 1) / (1 - a^J" 1 (*)^1 -1 (*)**■ ( 138 ) 

The integral in (138) may be evaluated by using (136) and the recur- 
rence relation (135), with I replaced by n also. Then, from (60), it is 
found that 

„ m _ n(n + l)(m 2 + rc 2 + n- 1) 
a,d (2n-l)(2n + 3) '" 

(n + 2)(n + 3)(n + \m\ + l)(n + \m\ + 2) 

2(2n + 3)(2n + 5) ' ,n+2 

(» - 2)(n - l)(n - lm|)(n - \m\ - 1) n , q v 
2(2ft-3)(2n-l) 5l - 2 - (139) 

Finally, from (61) and (133), 

d " 1 = \ f\W ~ ^ ^dx^ dPl dl (X) + m2p l m ^ x ) p " ml w] dx - ( 14 °) 
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The integral in (140) may be evaluated by using (136) and the relation 14 

(2ft + 1)(1 _^^w 

= (k + l)(/c + \m\)Pll\(x) - k(k - \m\ + l)P\$i(x). (141) 
Then from (61) it is found that 



m 2 n[(n -f- l) 2 — m 2 2 



n(n + l) ' (n+ l)(2n + l)(2n + 3) 



i«j - g < 

(w + 1) (n 2 - m 2 ) 1 _ (n + 3)(w + [m| + l)(n + |ro| + 2) 
+ >i(4n 2 - 1) ( l " 2(n + 1) (2n + 3) (2n + 5) 

v a (w-2)(n- lm|)(n- \m\ - 1) ,-.„. 

X Sl - +2 2n(2n-3)(2n-l) ai — *• (142) 

We note that the above results are consistent with the expressions 
given by Oguchi, 4 without derivation, for the integrals in (131), (61), 
and (62), subject to (133). 
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